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ABSTRACT 

We have identified outflows and bubbles in the Taurus molecular cloud based on the ~ 100 
deg 2 Five College Radio Astronomy Observatory 12 CO(1-0) and 13 CO(1-0) maps and the Spitzer 
young stellar object catalogs. In the main 44 deg 2 area of Taurus we found 55 outflows, of which 
31 were previously unknown. We also found 37 bubbles in the entire 100 deg 2 area of Taurus, all of 
which had not been found before. The total kinetic energy of the identified outflows is estimated 
to be ~ 3.9 x 10 45 erg, which is 1% of the cloud turbulent energy. The total kinetic energy of 
the detected bubbles is estimated to be ~ 9.2 x 10 46 erg, which is 29% of the turbulent energy 
of Taurus. The energy injection rate from outflows is ~ 1.3 x 10 33 erg s , 0.4 - 2 times the 
dissipation rate of the cloud turbulence. The energy injection rate from bubbles is ~ 6.4 x 10 33 
erg s _1 , 2-10 times the turbulent dissipation rate of the cloud. The gravitational binding 
energy of the cloud is ~ 1.5 x 10 48 erg, 385 and 16 times the energy of outflows and bubbles, 
respectively. We conclude that neither outflows nor bubbles can provide enough energy to 
balance the overall gravitational binding energy and the turbulent energy of Taurus. 
However, in the current epoch, stellar feedback is sufficient to maintain the observed turbulence 
in Taurus. 

Subject headings: ISM: jets and outflows - ISM: bubbles - Physical Data and Processes: turbulence - 
ISM: individual objects(Taurus) - ISM: kinematics and dynamics - Astronomical Databases: surveys 


1. Introduction 


Stars during their early stage of evolution ex¬ 
perience a ph ase of mass loss driven by strong 
stellar winds (Lada] ll985h . The stellar winds 
can entrain and accelerate ambient gas and in¬ 
ject momentum and energy into the surrounding 
environment, thereby significantly affect the dy¬ 
namics and structure of their parent molecular 


clouds (|Naravanan et al.l l2008t lArce et al.1 1201111 . 
Both outflows and bubbles are manifestations of 
strong stellar winds dispersing the surrounding 
gas. In general, collimated jet-like winds from 
young embedded protostars usually drive powerful 
collimated outflows, while wide-angle or spherical 
winds from the pre-main-sequence stars are more 
likely to drive less-collimated outflows or bubbles 
(lArce et al. 2011). A bubble is a partially or fully 
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enclosed three-dimensional stru cture whose pro¬ 


jectio n is a partial or full ring (jChurchwell et al 
2006). 


The kinetic energy of an outflow is very large 
(10 43 -10 48 erg; Ladal 1985 : Bachiller 1996h . im¬ 
plying a substantial input of mecha nical energy 


into its parent molecular cloud ( Solomon et al 
19811) . Feedback from young stars has been 


proposed as a signi ficant aspect of self-regulation 
of st ar formation ( Norman fe S ilk 198(J Francol 
1983h . Feedback may maintain the observed 
turbulence in molecular clouds and it may also 
be responsible for sta bilizing th e clouds against 
gravitational collapse ( Shu et ahlll987 ). The im¬ 
pact of outflows on surrounding gas has been 
studi ed p rimarily in small regions such _as Orion 
KL (Kwan fe Scoville 1976h L1551 (Snell et al 


1980h and GL 490 ( Lada fe Harvevl[l98ll) on scales 

less than 10'. Recently there have been a few 
studies rela ted to outflow feedback in nearby 
clouds. Arce et all ( 2010h undertook a complete 
survey of outflows in Perseus and found that 
outflows have an important impact on the environ¬ 
ment immediately surrounding localized regions 
of active star formation, but that outflows have 
insufficient energy to feed the obs erved turbulence 


i n the entire Perseus complex. iNakamura et al 
( 2011a ) and Nakamura et al. ( 2011bl) studied the 


outflows in the p Ophiuchi main cloud and Ser¬ 
pens south, respectively. Both studies concluded 
that outflows can power the supersonic turbu¬ 
lence in their parent molecular cloud but do not 
have enough momentum to support the entire 
cloud against the global gravitational contraction. 
Narayanan et al. ( 2012 ) identified 20 outflows in 
the Taurus region and concluded that outflows 
cannot sustain the observed turbulence seen in the 
entire cloud. In this paper, we report a systematic 
and detailed search for outflows around sources 
from the Spitzer Space Telescope (hereinafter 
referred to as Spitzer) young stellar object (YSO) 
catalog and then estimated their impact on the 
entire Taurus molecular cloud. 

Similar to outflows, bubbles are important 
morphological features in star formation pro¬ 
cess, which can give information about spher¬ 
ical stellar winds and physical properties of 


their surrounding environments (jChurchwell et al 


2006) . Parsec-scale bubbles ar e us ually found in 
massive star-forming regions (IHever et al.1 1_1992; 


Churchwell et al. 2006 1 

2007; 

Beaumont & Williams 

2010) Dcharveng et al. 

2010 

). The conventional 


thought has been that high-mass stars can drive 
spherical winds and easily create the observed 
bubbles, while the spherical winds from low- and 
intermediate-mass sta rs are too weak to produce 
bubbles. However, Arce et ahl ( 2011 ) studied 


shells (bubbles) in Perseus, a nearby low-mass 
star-forming molecular cloud, and concluded that 
the total energy input from outflows and shells is 
sufficient to maintain the turbulence. 

The Ta urus molecular clo ud is at a distance 
of 140 pc ( Torres et al • 2001 . It_ covers an area 
of m ore than 100 deg 2 (jlJngerechts fc Thaddeus 
1987 ). Using the J=2-l line of 12 CO, 13 
outflows have been found around l ow-mass 
embe dded YSOs in Taurus dBontemps et al. 
1996 ). There are 13 high velocity molecular 


outfl ows in Taurus included in the catalog of 


Wu e t al. (120041) . Using JCMT-HARP 12 CO 
J=3-2 observations, 16 outflows have been found 
L1495, a ‘bowl-sha ped’ region in the NW 


corner of Taurus ( Davis et al. 2010h . Recently, 20 
outflows have been identified, 8 of which were new 
detections with the Five College Radio Astronomy 
Observatory (FCRAO) 12 CO J=l-0 and 13 CO 
J=l-0 data cu bes covering the entir e Taurus 
molecular cloud (jNaravanan_et_ah 2012 ). Th e up- 
to-date catalog of YSOs dRebull et al.ll2010i ) from 
the Spitzer provides an opportunity to search for 
outflows and bubbles in a more comprehensive 
manner. Here we present a systematic and 
detailed search for outflows and bubbles in the 
vicinity of young stellar objects and estimate their 
impact on the overall Taurus molecular cloud. 

The paper is organized as follows. In § [2] 
we describe the data used in the study. The 
details including searching methods, morphology 
and physical parameters of outflows and bubbles 
are presented in § OH and § HI respectively. The 
driving sources of outflows and bubbles, their 
energy feedback to the parent cloud and the 
related comparison between Taurus and Perseus 
are discussed § 0 In § [G] we summarize the main 
results. 

2. The Data 

In our study we used the 12 CO(1-0) and 
13 CO(1-0) data observed with 13.7 m FCRAO 
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telescope ( Narayanan et al.l 12008 ). We also 
adopted the up-to-date catalog of Spitzer YSOs, 
where 215 YSOs and 140 new YSO can didates in 
Taurus are reported ( Rebull et al.l 2010h . 


2.1. FCRAO CO Maps 

The FCRAO CO survey was taken between 
2003 and 2005. The 12 CO and 13 CO maps 
are centered at a=04 h 32 m 44.6 s , 6=24°25'44.6" 
(J2000) covering an area of approximately 100 
deg 2 . The full width at half maximum (FWHM) 
beam width is 45" for 12 CO and is 47" for 13 CO. 
The pixel size of the resampled data is 20", which 
corresponds to 0.014 pc at a distance of 140 pc. 
There are 80 channels for 12 CO and 76 channels 
for 13 CO, covering approximately —5 to +14.9 km 
s _1 . The width of a velocity channel is 0.254 
km s _1 for 12 CO and 0.266 km s _1 for 13 CO 
( Narayanan et al. 20081 : [Goldsmith et al. 2008 ). 


2.2. Spitzer MIPS Images 

The MIPS (Multi-band Imaging Photometer 
for Spitzer; Rieke et al. 20041) maps were created 
as part of the final products from t he Spitzer 


Legac y Taurus I and II surveys (jPadgett et al 
20071 ). The data were obtained in fast scan 
mode in three bands: 24, 70, and 160 /im, 
over an area of 44 deg 2 . The observations were 
performed in three epochs between 2005 and 
2007, with an integration time of 30 s (24 /im) 
and 15 s (70 & 160 /im). The maps were 
created using the basic calibrated data (BCDs) 
and coadded using the Spitzer software package 
MOPEX (Mosaick ing and Point Source Extractor; 
Makovoz fc Marleaul 120051) . Despite the fact that 
the data were taken with interleaved scan legs to 
provide optimal coverage at 70 and 160 /rm, some 
small gaps remained, in particular at 160 /mi. To 
mitigate this effect, the 160 /mi final mosaic was 
created using 32 arcsec pixels, instead of the native 
16 arcsec/pixel scale. This pixel scale matches 
quite well with the ~40 arcsec beam at 160 /(m 
wavelength. The 24 and 70 /mi maps adopted the 
standand 2.5 and 4 arcsec/pixel scale, respectively, 
to properly sample their respective 6 and 18 
arcsec beams. The maps were used successfully for 
photometric purposes to ide ntify new sources i n 
the Taurus Molecular Cloud ( Rebull et al • 2010). 


3. Outflows 

We identified 55 outflows around the Spitzer 
YSOs in the 44 deg 2 area of Taurus. In total 31 of 
the detected outflows were previously unknown. 
In the following subsections, we describe the 
searching procedure of outflows, the morphology, 
physical properties, and the comparison between 
our findings and the known ones. 

3.1. The Search Procedures for Outflows 

Instead of a blind search, we focused on seeking 
outflows around YSOs. The search procedure 
was performed with an Interactive Data Language 
(IDL) pipeline. We plotted spectra, position 
velocity diagrams (hereafter P-V diagrams) and 
integrated intensity maps to identify the outflows 
around the 355 YSOs which Spitzer identified 
in Taurus. Detailed steps of the search are the 
following. 

(1) We plotted 12 CO contours (hereafter con¬ 
tour map) overlaid on a 13 CO grey-scale image 
around a YSO. According to the scale and velocity 
range of previously-detected outflows in Taurus, 
we chose two sizes (10' and 20') and three sets of 
velocity intervals (-1 to 3.5 km s 1 , -1 to 4.5 km 
s _1 and -1 to 5.5 km s _1 for blue; 7.5 to 13 km 
s _1 , 8.5 to 13 km s _1 and 9.5 to 13 km s _1 for red) 
to plot the contour maps. We plotted the maps 
with 3 sets of velocity intervals and two scales 
around the 355 YSOs automatically. In total, 2130 
maps were obtained. We inspected these maps 
to identify outflow candidates according to the 
morphology of the blue and red lobes. In the end, 
74 candidates were selected. 

(2) We plotted 12 CO P-V diagrams along four 
directions (at position angles0 of 0°, 45°, 90° and 
135°) on three scales (20', 40' and 60') around the 
74 candidates. The size and high-velocity range 
of outflow candidates were determined roughly by 
checking the P-V diagrams. When the velocity 
bulge appears along the direction away from the 
central velocity, we marked it as the start of 
the high-velocity wing. And along the above 
direction, the maximum velocity corresponding to 
the outermost contour is the end of this high- 


1 The position angle is defined as the angle measured from 
the north clockwise to the direction along which we plotted 
the P-V diagram. 
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velocity wing. We further confirmed each of the 
diagrams individually by visual inspection. The 
position range of the entire high-velocity bulge 
along the position axis was considered to be the 
rough size of the outflow. If more than one central 
velocity is found in the P-V diagram, it likely 
has multiple velocity components ( Wu et ahll2005 ) 
and thus will be excluded from the list of outflow 
candidates. Therefore, 19 candidates with multi¬ 
ple velocity components were eliminated and the 
remaining 55 outflow candidates were considered 
to be possible outflows. 

(3) Using the rough sizes and velocity ranges 
obtained in step (2), we plotted contour maps for 
the remaining 55 outflows. P-V diagrams were 
plotted through the midpoint of the blue and red 
peaks (bipolar outflow) or through the peak of the 
lobe in the case of monopolar outflow, at position 
angles spaced by 15°. We chose the angle with 
the most prominent bulge along the velocity axis 
to determine the velocity interval of the outflows. 
Then we plotted the contour map again with this 
velocity interval. 

(4) Finally we plotted the average spectra of the 
blue and red lobes. According to the morphology 
of P-V diagrams and contour maps, we divided 
the outflows into five classes. The higher is 
the ranking, the more likely it is that we have 
identified an outflow. We define a typical P-V 
diagram (TPV) and a representative contour map 
(RCM) as follows. If there is obvious high velocity 
gas which can be seen by the protuberance along 
the velocity axis on the P-V diagram and the 
high velocity range is not less than 1 km s , 
then we regard the P-V diagram as a TPV. If the 
outermost contour of the lobe is closed and we 
can see a clear and unbroken lobe on the contour 
map, then we regard the contour map as a RCM. 
Table [T| shows our criteria for outflow classification 
with an “x” meaning that it satisfies a certain 
condition. Having both 12 CO TPV and 12 CO 
RCM is required for high ranking (Class A + ), 
but having 13 CO TPV or 13 CO RCM gives a 
lower ranking (Class A - and Class B~) because 
13 CO is generally optically thin in outflows and 
we usually found lobes of outflows with 12 CO 
not 13 CO. Having only 12 CO RCM was divided 
into the lowest ranking (Class C + ) because the 
high velocity gas in P-V diagram is not obvious. 
The primary condition to identify an outflow is 


having high velocity gas which can be seen from 
the protuberance along the velocity axis on the 
P-V diagram. 

3.2. The Results of Outflow Search 

Following the steps in § 13.1 1 we found 55 outflows 
in the 44 deg 2 area of Taurus molecular cloud. All 
the outflows we detected were listed in Table [2] 
Each outflow is referred as a “Taurus Molecular 
Outflow” (TMO). We present the locations, po¬ 
larities and scales of the outflows overlaid on the 
Spitzer MIPS Image in Fig. [Q There are 31 new 
ones among all the detected outflows. We have 
thus increased the total number of known outflows 
by a factor of 1.3. 


Table 1: Criteria for Outflow Classification and 
Outflow Distribution 


Class 

i2 CO 

TPV 

i2 co 

RCM 

i3 CO 

TPV 

i3 co 

RCM 

Outflow 

Numbers 

Percentage 

A+ 

X 

X 



24 

43.6% 

A" 

X 

X 

X 

X 

18 

32.7% 

B+ 

X 




4 

7.3% 

B~ 

X 


X 


1 

1.8% 

C+ 


X 



8 

14.5% 


4 







Table 2 

Outflows in Taurus 


Outflow 

Name 

R.A. 

(J2000) 

Dec. 

(J2000) 

Common 

Name 

YSO 

Type a 

Outflow 

Class 

Po. b 

New 

Detection 0 



Ref. 

TMO.Ol 

04 

11 

59.7 

29 

42 

36 

_ 

III 

A+ 

MR 

N 

1 



TMO.02 

04 

14 

12.2 

28 

08 

37 

IRAS 04113+2758 (L1495) 

I 

A+ 

Bi 

N 

1, 

2, 

3 

TMO.03 

04 

14 

14.5 

28 

27 

58 

- 

II 

A+ 

Bi 

Y 

4 



TMO.04 

04 

18 

32.0 

28 

31 

15 

- 

flat 

A + 

Bi 

Y 

4 



TMO.05 

04 

19 

41.4 

27 

16 

07 

IRAS 04166+2706 

I 

A+ 

MB 

N 

1, 

2, 

5, 6, 7 

TMO.06 

04 

19 

58.4 

27 

09 

57 

IRAS 04169+2702 

I 

A+ 

Bi 

N 

1, 

2, 

3, 7 

TMO.07 

04 

21 

07.9 

27 

02 

20 

IRAS 04181+2655 

I 

A+ 

Bi 

N 

2, 

3, 

6, 7 

TMO.08 

04 

22 

15.6 

26 

57 

06 

FS Tau B 

I 

A+ 

Bi 

N 

1, 

2 


TMO.09 

04 

23 

25.9 

25 

03 

54 

- 

II 

A + 

MR 

Y 

4 



TMO.IO 

04 

24 

20.9 

26 

30 

51 

- 

II 

A + 

Bi 

Y 

4 



TMO.ll 

04 

24 

45.0 

27 

01 

44 

- 

III 

A+ 

MR 

Y 

4 



TMO.12 

04 

29 

30.0 

24 

39 

55 

Haro 6-10 

I 

A+ 

MR 

N 

1, 

6, 

8 

TMO.13 

04 

31 

10.4 

25 

41 

29 

- 

- 

A+ 

MR 

Y 

4 



TMO.14 

04 

31 

58.4 

25 

43 

29 

- 

Ill 

A+ 

Bi 

Y 

4 



TMO.15 

04 

32 

14.6 

22 

37 

42 

- 

flat 

A+ 

Bi 

Y 

4 



TMO.16 

04 

32 

31.7 

24 

20 

02 

L1529 

II 

A + 

Bi 

N 

1, 

6, 

9, 10 

TMO.17 

04 

32 

32.0 

22 

57 

26 

IRAS 04295+2251 (L1536) 

I 

A + 

MR 

N 

3, 

7 


TMO.18 

04 

32 

43.0 

25 

52 

31 

- 

II 

A+ 

Bi 

Y 

4 



TMO.IO 

04 

34 

15.2 

22 

50 

30 

- 

II 

A+ 

MR 

Y 

4 



TMO.20 

04 

37 

24.8 

27 

09 

19 

- 

III 

A+ 

MR 

Y 

4 



TMO-21 

04 

39 

53.9 

26 

03 

09 

L1527 

I 

A+ 

Bi 

N 

1, 

3, 

6, 7, 11, 12, 13 

TMO_22 

04 

41 

08.2 

25 

56 

07 

IRAS 04381+2540 (TMC-1) 

flat 

A + 

MB 

N 

1, 

6, 

7, 14 

TMO-23 

04 

41 

12.6 

25 

46 

35 

IRAS 04381+2540 (TMC-1) 

I 

A+ 

MR 

N 

1, 

6, 

7, 14 

TMO-24 

04 

42 

07.7 

25 

23 

11 

IRAS 04390+2517 (LkHa 332) 

II 

A+ 

MR 

N 

3 



TMO.25 

04 

18 

58.1 

28 

12 

23 

IRAS 04158+2805 (L1495) 

flat 

A - 

MB 

Y 

4, 

7 


TMO.26 

04 

23 

18.2 

26 

41 

15 

- 

II 

A” 

MR 

Y 

4 



TMO-27 

04 

26 

56.2 

24 

43 

35 

IRAS 04239+2436 (HH 300) 

I 

A~ 

MR 

N 

1, 

3, 

6, 15 

TMO_28 

04 

27 

02.6 

26 

05 

30 

IRAS 04240+2559 (DG Tau) 

I 

A" 

Bi 

N 

1, 

16 

TMO-29 

04 

27 

02.8 

25 

42 

22 

- 

II 

A“ 

Bi 

Y 

4 



TMO.30 

04 

27 

57.3 

26 

19 

18 

IRAS 04248+2612 

flat 

A" 

MR 

N 

1, 

3 


TMO_31 

04 

28 

10.4 

24 

35 

53 

- 

flat 

A - 

MR 

Y 

4 



TMO.32 

04 

30 

51.7 

24 

41 

47 

IRAS 04278+2435 (ZZ Tau IRS) 

flat 

A” 

MR 

N 

1, 

6, 

17 

TMO.33 

04 

32 

15.4 

24 

28 

59 

IRAS 04292+2422 (Haro 6-13) 

flat 

A’ 

Bi 

N 

1, 

3 


TMO_34 

04 

33 

07.8 

26 

16 

06 

- 

III 

A - 

MB 

Y 

4 



TMO_35 

04 

33 

10.0 

24 

33 

43 

- 

III 

A“ 

MR 

Y 

4 



TMO.36 

04 

33 

16.5 

22 

53 

20 

IRAS 04302+2247 

I 

A“ 

Bi 

N 

1, 

3, 

7 

TMO_37 

04 

33 

34.0 

24 

21 

17 

- 

II 

A - 

MR 

Y 

4 



TMO.38 

04 

33 

36.7 

26 

09 

49 

- 

II 

A~ 

MB 

Y 

4 



TMO_39 

04 

35 

57.6 

22 

53 

57 

IRAS 04328+2248 (HP Tau) 

II 

A’ 

Bi 

N 

3 



TMO.40 

04 

39 

11.2 

25 

27 

10 

HH706 

- 

A” 

MR 

N 

1 



TMO.41 

04 

39 

13.8 

25 

53 

20 

IRAS 04361+2547 (TMR-1) 

I 

A“ 

Bi 

N 

1, 

3, 

6, 7, 18 

TMO_42 

04 

48 

02.3 

25 

33 

59 

Tau A 8 

III 

A“ 

Bi 

Y 

4 



TMO_43 

04 

18 

51.4 

28 

20 

26 

HH156 

I 

B+ 

MB 

Y 

2 



TMO.44 

04 

20 

21.4 

28 

13 

49 

- 

flat 

B+ 

Bi 

Y 

4 



TMO.45 

04 

26 

53.3 

25 

58 

58 

- 

I 

B+ 

Bi 

Y 

4 



TMO.46 

04 

39 

56.1 

26 

28 

02 

- 

- 

B+ 

Bi 

Y 

4 



TMO.47 

04 

35 

35.3 

24 

08 

19 

IRAS 04325+2402 (L1535) 

I 

B“ 

Bi 

N 

1, 

3, 

6, 17, 19, 20 

TMO.48 

04 

15 

35.6 

28 

47 

41 

- 

I 

c+ 

Bi 

Y 

4 



TMO_49 

04 

17 

33.7 

28 

20 

46 

- 

II 

C+ 

MR 

Y 

4 



TMO.50 

04 

18 

10.5 

28 

44 

47 

- 

I 

C+ 

MR 

Y 

4 



TMO.51 

04 

18 

31.1 

28 

16 

29 

- 

II 

c+ 

MR 

Y 

4 



TMO.52 

04 

18 

31.2 

28 

26 

17 

- 

I 

c+ 

MR 

Y 

4 



TMO.53 

04 

18 

41.3 

28 

27 

25 

- 

flat 

c+ 

MR 

Y 

4 



TMO.54 

04 

21 

54.5 

26 

52 

31 

- 

II 

c+ 

Bi 

Y 

4 



TMO.55 

04 

29 

04.9 

26 

49 

07 

IRAS 04260+2642 

I 

c+ 

MR 

N 

4 




a The YSO classification from lRebull et al.1 (12010 1. “flat” represents a flat-spectrum YSO, which is intermediate between Class I and II. 
b The polarity of the outflow. “Bi” represents bipolar outflow, “MB” and “MR” indicate blue and red monopolar outflow, respectively. 
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The column represents whether the outflow is detected for the first time in our study. Y = new, N = has reported in previous work. 


References. — (1) Narayanan et al. (2012); (2) Davis et al. (2010); (3) Moriarty-Schieven et al. (1992); (4) this paper; (5) Tafalla et al. 
(2004); (6) Wu et al. (2004); (7) Bontemps et al. (1996); (8) Stojimirovic, Narayanan 8z Snell (2007); (9) Lichten (1982); (10) Goldsmith et 
al. (1984); (11) Tamura et al. (1996); (12) Hogerheijde et al. (1998); (13) Zhou, Evans &; Wang (1996); (14) Chandler et al. (1996); (15) 
Arce Sz Goodman (2001); (16) Mitchell et al. (1994); (17) Heyer, Snell Sz Goldsmith (1987); (18) Terebey et al. (1990); (19) Myers et al. 
(1988); (20) Wu, Zhou, &; Evans II (1992). 
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Table m lists the numbers and percentages in 
the five classes of outflows. We can see Class A + 
and class A - account for 76.3% of all the detected 
outflows. These two types can be considered as 
the “most probable” outflows in our study. Table 
[3] lists the numbers of previously known and newly 
detected outflows in different classes. We found 
more new outflows of Class A + and class A - , 
which account for 64.5% of all the newly detected 
outflows. That is, most of the new outflows we 
found are likely true outflows. 

Table m lists the outflow numbers and per¬ 
centages according to the types of their driving 
sources. Class I accounts for 36.4%, which is 
the largest proportion of all the YSOs driving 
outflow. The outflows driven by the Class I YSOs 
are closer to the YSOs and have more collimated 
bi-polar morphology. Compared with the Class I, 
Class III YSOs drive a small proportion (12.7%) 
of outflows, which tend to be farther from the 
YSOs. This indicates that the outflows from 
Class III YSOs are more evolved than those from 
Class I YSOs. We also found three outflows 


(TMO_13, TMO_40, TMO_46) without YSOs, 
indicating that they are possibly Class 0 objects. 
Among the three outflows, TMO_13 and TMO_40 
are newly found in our study, whi l e TM O_46 has 
been reported in Narayanan et ajJ (12012 '). 


3.3. Morphology of Outflows 


We found 25 bipolar, 22 monopolar redshifted 
and 6 monopolar blueshifted outflows. Bipo¬ 
lar and redshifted outflows account for the vast 
majority of outflows in the Taurus molecular 
cloud. This is consis tent with the results of 
Narayanan et al.1 ( 20121 ). Fig. [5] - Fig. [59] show 
the 12 CO integrated intensity map, 12 CO P-V 


diagram and average spectrum for each outflow. 
For Class A - and B~ outflows, we also plotted 
the 13 CO integrated intensity maps and 13 CO P- 
V diagrams. 


3.4. Comparison with Previously Found 
Outflows 


Using the FCRAO large-scale survey data 


(Narayanan et al.1 2008) _and the latest YSO 


catalog from Rebull et al.1 ( 2010h we were able 
to identify the previously known outflows, obtain 
more complete morphology, and find additional 
new outflows. The YSO catalog is also convenient 
for identifying the driving sources of the outflows. 
Comparing with the previous works, we confirmed 
more driving sources of outflows. 

L1527 (TMO_21) is a typical outflow in Taurus 
( Narayanan et al. 2012 ; Hogerheiide et al. 19981 ; 


Zhou et al.l ll996l) . The P-V diagram and contour 


map in ou r wo rk are v ery simil ar t o th ose in 
Hogerheiide et al. ( 1998 1 and Narayanan et al 


(hom TMC L08 fSST 042215.6+265706) and FS 


Tau B in Narayanan et al. ( 20121 ) are the same 
outflow with the same location. They have the 
same structure, which can be seen in o ur Fig . 
m and Figure 15 in iNaravanan et al.1 (120121) . 

In addition, TMCL30 (SST 042757.3+261918), 
TMCL32 (SST 043051.7+244147), TMCL33 (SST 
043215.4+242859) and TMCL41 (SST 043913.8+255320) 
also have the same morphology as IRAS 04248+2612, 

ZZ Tau IRS, IRAS 04292+2 4 22 an d IRAS 
04361+2547 in INaravanan et al.1 ( 2012 ). respec¬ 
tively. These confirm the general consistency 
between the two works in terms of strong and 
extended outflows. 

For TMCL02 (SST 041412.2+280837), we ob¬ 
tained a good bipolar structure shown in the upper 


Table 3: The Numbers of Outflows in Different 
Classes _ 


Type 

Previously 

Known 

Newly 

Detected 

A+ 

13 

11 

A- 

9 

9 

B+ 

0 

4 

B“ 

1 

0 

C+ 

1 

7 


Table 4: Number 
Types of YSOs 

of Outflows 

Around Different 

YSO 

Type 

Number of 
Outflows 

Percentage 

I 

20 

36.4% 

Flat 

10 

18.2% 

II 

15 

27.3% 

III 

7 

12.7% 

No YSO 

3 

5.4% 


7 




















































































30° 


29° 


28° 


27° 

o 

O 26° 

<N 

Z 25° 

CD 

Q 

24° 


23° 


22 ° 


21 ° 

04h32m 


24m 20m 16m 

RA (J2000) 



Fig. 1. — Outflows and bubbles overlaid on the Spitzer 160 fin 1 (grey), 70 fim (green) and 24 /mi (red) image 
from the Spitzer Legacy Taurus I and II surveys. The blue and red sectors represent blue and red lobes of 
outflows in Table [U respectively. The radius of the sectors shows the scale of the outflow, while the direction 
of the sectors shows the projected direction of the outflow. The green rings and arcs represent the expanding 
and broken bubbles, respectively. The thickness and radius of the arcs and rings are the actual thickness 
and radius of the bubble structures. T he ope n orange circles show the locations of YSOs, which were listed 
in Table 6 and Table 7 of lRebull et al. ( 2010h . 









left panel of Fig. [5J while Narayanan et al.i ( 20081) 


considered this outfl ow (IRAS 0 4113+ 2758) to 
be redshifted only. iDavis et~ahl ( 20 id ) did not 


identify the driving source of this outflow (named 
by W-CO-flowl), while we determined that the 
YSO SST 04141 2 .2+280 837 is driving the outflow. 
Moriartv-Schieven et al. ( 1992 1 only presented the 


central spectrum of IRAS 04390+2517 and IRAS 
04328+2248, while we illustrated the two outflows 
(TMO_24 and TMO_39) more clearly through 
contour maps and P-V diagrams. 

The morphology of TMCL07 (SST 042107.9+270220) 


of their small sizes. 

3.5. Physical Parameters of Outflows 

To study the effects of outflows on their en¬ 
vironment we calculated their masses, momenta, 
kinetic energy and energy deposition rates. The 
total column density of the outflowing gas is 


AtotrCO) = 


3 k 2 T e . 


47t 3 fJ^hv 2 exp (—hv/kTex) 


T s dv, 


(1) 

where k = 1.38 x 10 16 erg/K, h = 6.626 x 10 27 
erg s, lia = 0.112 x 10 -18 esu, v = 115.2712 x 


in Pavis et al.l 

(2010). This outflow 

was also 

10 9 Hz and T s is the observed source antenna 

reported by 

Moriartv-Schieven et al. 

( 

1992 

), 

)■ 

temperature with proper correction for antenna 

Bontemps et al. 

199C}) and Wu et al 

( 

2004 

efficiency. We assumed an excitation temperature 

However. Naravanan et al. ( 

2012) did not find it 

of 25 K. The excitation temperature assumed in 


with the same FCRAO survey data. iMitchell et al 


(119941 ) reported that IRAS 04240+2559 was a 
monopolar redshifted outflow with 12 CO (3-2) 
line. But at this location we found the well-defined 
bipolar outflow as shown in Fig. [32] 

As for L1529, Lichten ( 19821 ) presented high- 


velocity 12 CO wings observed by antenna No. 2 of 
the C altech 10.4 m array, but iGoldsmith et al.l 
( 19841) did not find any high-velocity gas in 


observations at FCRAO. We identified a bipolar 
outflow named TMO_16 (SST 04 3231.7+ 242002) 
and demonstrated the result of Lichtenl (.1 982 ) 



found no outflow. We found a red monopolar 
outflow as shown in Fig. l2ll 

We have found 31 new outflows which are 
labeled “Y” in the eighth column of Table [2] 
Two of these new outflows wer e not identified 
as outflows in the literature. Bontemns et al.l 


(119961) considered IRAS 04158+2805 (L1495) but 


did not find any sign of outflow activity in the 
12 CO (2-1) transition at the location of TMO_025 
(SST 041858.1+281223). IPavis et al.1 (l2010h had 


some doubt about CO flow of CoKU Tau-1 when 
analyzing the 12 CO (3-2) emission, while we found 
TMCL043 (SST 041851.4+282026) at this site. 
The rest of the new outflows have not been 
reported in the literature and are identified as 
outflows for the first time. All of the new outflows 
are of small angular extent, less than 10'. They 
were missed in previous searches perhaps because 


the literature (Zhouetal. 

1996 

Tarnura et al. 

1996; Ohashi et al. 1997a 

d: Pavis et al. 2010; 

Naravanan et al. 2012) ranges from 11 K to 50 


K. The lowest temperature will decrease the 
mass estimate by a factor of 3 and the highest 
temperature will increase the mass estimate by a 
factor of 2.2. The detailed derivations regarding 
the physical parameters of the outflows are given 
in appendix A. 


Arce & Goodman ( 20011 ) described three major 


issues that can cause uncertainties in the calcula¬ 
tion of an outflow’s parameters, namely, the incli¬ 
nation, opacity, and blending. Our prescriptions 
are the following. (1) We defined the inclination 
angle of the outflow as the angle between the 
long axis of the outflow and the line of sight. 
Since the outflows with small inclination angle 
(especially when the outflow is perpendicular 
to the plane of the sky) are hard to detect, our 
outflow searching is biased to those with large 
inclination angle. If the inclination angle 9 is 
randomly distributed, the average value is given 

by 

r /2 

(9)= 0sin 9d6. (2) 

Jo 

From the above formula we got the average in¬ 
clination angle of 57.3°, which differs from the 
usually used median value of 45°. Then the 
velocity and the dynamic age, fdyn, should be 
scaled up by a factor of 1.9 and 0.64, respectively. 
(2) Using the 12 CO and 13 CO data we can correct 
for the opacity in the 12 CO line when 12 CO 
emission of an outflow is optically thick. The 


9 














































































































algorithm for the opacity correction is described 
in appendix A. (3) We probably missed some 
low-velocity outflowing gas, which blended into 
the ambient gas, when we conservatively deter¬ 
mined t he emission only from outflows. Previous 
studies (iMareulis fc Lada Il985t lArce et al1l2010t 


Narayanan et al.1 120121 ) showed that neglecting 


this gas results in the underestimate of the outflow 
mass almost by a factor of 2. 

Table [5] gives the length, mass, momentum, 
kinetic energy, dynamical timescale and the lumi¬ 
nosity of the outflows in Taurus. 
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Table 5 

Physical Parameters of Outflows 


Outflow 

Name 

Lobe 

^avg 

(km s _1 ) 

Area b 

(arcmin) 

Length 

(pc) 

Mass 

(Mq) 

Momentum 
( Mq km s —1 ) 

Energy 
(10 43 erg) 

^•dyn 

(10 5 yr) 

^flow 

(10 30 erg s -1 ) 

TMO.Ol 

Blueshifted 

_ 

_ 

_ 

_ 

_ 

_ 

_ 

_ 


Redshifted 

2.5 

11 x 25 

1.11 

0.083 

0.205 

0.51 

4.4 

0.37 

TMO.02 

Blueshifted 

3.4 

8 X 12 

0.58 

0.168 

0.570 

1.92 

1.7 

3.63 


Redshifted 

3.6 

6x5 

0.31 

0.061 

0.220 

0.79 

0.8 

2.97 

TMO.03 

Blueshifted 

3.1 

6x9 

0.41 

0.050 

0.156 

0.49 

1.3 

1.19 


Redshifted 

2.4 

11 X 2 

0.46 

0.013 

0.030 

0.07 

1.9 

0.12 

TMO.04 

Blueshifted 

2.6 

6x8 

0.41 

0.026 

0.069 

0.18 

1.5 

0.36 


Redshifted 

2.5 

3x8 

0.33 

0.005 

0.013 

0.03 

1.3 

0.08 

TMO.05 

Blueshifted 

4.6 

9 x 18 

0.83 

0.078 

0.355 

1.61 

1.8 

2.89 


Redshifted 

- 

- 

- 

- 

- 

- 

- 

- 

TMO.06 

Blueshifted 

2.8 

5x5 

0.29 

0.038 

0.109 

0.31 

1.0 

0.96 


Redshifted 

3.9 

3x5 

0.21 

0.010 

0.038 

0.15 

0.5 

0.88 

TMO-07 

Blueshifted 

2.1 

7x7 

0.38 

0.012 

0.025 

0.05 

1.8 

0.09 


Redshifted 

2.7 

2x2 

0.13 

0.002 

0.005 

0.01 

0.5 

0.08 

TMO.08 

Blueshifted 

3.1 

5x4 

0.26 

0.013 

0.039 

0.12 

0.8 

0.48 


Redshifted 

2.8 

2x2 

0.12 

0.001 

0.003 

0.01 

0.4 

0.07 

TMO.09 

Blueshifted 

- 

- 

- 

- 

- 

- 

- 

- 


Redshifted 

1.8 

17 x 18 

1.01 

0.311 

0.565 

1.02 

5.5 

0.59 

TMO.10 

Blueshifted 

2.3 

3x7 

0.29 

0.010 

0.023 

0.05 

1.3 

0.13 


Redshifted 

2.9 

7x3 

0.30 

0.016 

0.047 

0.14 

1.0 

0.43 

TMO.ll 

Blueshifted 

- 

- 

- 

- 

- 

- 

- 

- 


Redshifted 

2.6 

7x7 

0.40 

0.020 

0.052 

0.14 

1.5 

0.29 

TMO.12 

Blueshifted 

- 

- 

- 

- 

- 

- 

- 

- 


Redshifted 

2.5 

4x3 

0.22 

0.014 

0.036 

0.09 

0.8 

0.33 

TMO.13 

Blueshifted 

- 

- 

- 

- 

- 

- 

- 

- 


Redshifted 

2.2 

7x4 

0.34 

0.053 

0.115 

0.25 

1.5 

0.51 

TMO.14 

Blueshifted 

2.6 

5x5 

0.29 

0.019 

0.050 

0.13 

1.1 

0.38 


Redshifted 

3.2 

7x4 

0.31 

0.031 

0.099 

0.32 

0.9 

1.05 

TMO.15 

Blueshifted 

1.9 

3x3 

0.18 

0.013 

0.024 

0.05 

0.9 

0.16 


Redshifted 

1.5 

3x4 

0.20 

0.010 

0.016 

0.02 

1.3 

0.06 

TMO.16 

Blueshifted 

2.3 

2x2 

0.11 

0.005 

0.012 

0.03 

0.4 

0.21 


Redshifted 

2.0 

4x4 

0.23 

0.023 

0.046 

0.09 

1.1 

0.25 

TMO.17 

Blueshifted 

- 

- 

- 

- 

- 

- 

- 

- 


Redshifted 

1.7 

5 x 14 

0.60 

0.042 

0.072 

0.12 

3.5 

0.11 

TMO.18 

Blueshifted 

2.0 

10 x 3 

0.42 

0.012 

0.024 

0.05 

2.1 

0.07 


Redshifted 

1.9 

3x4 

0.20 

0.007 

0.013 

0.02 

1.0 

0.08 

TMO.19 

Blueshifted 

- 

- 

- 

- 

- 

- 

- 

- 


Redshifted 

2.2 

4x4 

0.23 

0.008 

0.017 

0.04 

1.0 

0.11 

TMO.20 

Blueshifted 

- 

- 

- 

- 

- 

- 

- 

- 


Redshifted 

4.0 

7x3 

0.31 

0.017 

0.067 

0.26 

0.8 

1.08 

TMO.21 

Blueshifted 

3.1 

4x2 

0.17 

0.004 

0.012 

0.04 

0.5 

0.22 


Redshifted 

2.9 

5x3 

0.22 

0.011 

0.032 

0.09 

0.7 

0.40 

TMO.22 

Blueshifted 

1.6 

10 x 22 

0.97 

0.226 

0.351 

0.54 

6.1 

0.28 


Redshifted 

- 

- 

- 

- 

- 

- 

- 

- 

TMO_23 

Blueshifted 

- 

- 

- 

- 

- 

- 

- 

- 


Redshifted 

3.7 

3x6 

0.27 

0.019 

0.072 

0.26 

0.7 

1.18 

TMO.24 

Blueshifted 

- 

- 

- 

- 

- 

- 

- 

- 


Redshifted 

3.2 

4x6 

0.29 

0.015 

0.047 

0.15 

0.9 

0.52 

TMO.25 

Blueshifted 

2.4 

5x5 

0.31 

0.267 

0.640 

1.52 

1.2 

3.87 


Redshifted 

- 

- 

- 

- 

- 

- 

- 

- 

TMO.26 

Blueshifted 

- 

- 

- 

- 

- 

- 

- 

- 


Redshifted 

2.4 

13 x 8 

0.60 

0.369 

0.880 

2.09 

2.5 

2.68 

TMO.27 

Blueshifted 

- 

- 

- 

- 

- 

- 

- 

- 


Redshifted 

2.4 

6x7 

0.37 

0.074 

0.174 

0.41 

1.5 

0.86 

TMO.28 

Blueshifted 

2.2 

4x6 

0.30 

0.170 

0.371 

0.80 

1.3 

1.91 


Redshifted 

3.9 

3x2 

0.15 

0.031 

0.122 

0.47 

0.4 

3.80 

TMO.29 

Blueshifted 

1.9 

9x4 

0.40 

0.264 

0.492 

0.91 

2.1 

1.37 


Redshifted 

2.1 

6x5 

0.32 

0.076 

0.161 

0.34 

1.5 

0.74 
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The distributions of length, mass, energy and 
dynamical timescale of outflows are shown in Fig. 
[2j The extents of outflows are in the range of 0.1 - 
1.11 pc. 79% of outflows are smaller than 0.6 pc. 
The mass of 54% of outflows is between 0.01 M 0 
and 0.1 Mg. The outflows with mass lower than 
0.01 M 0 and higher than 0.1 Mq account for 17% 
and 29% of the total, respectively. The energy of 
48% of outflows is in the range 10 42 - 10 43 erg. 
The outflows with energy lower than 10 42 erg and 
higher than 10 43 erg account for 31% and 21%, 
respectively. The dynamical timescales of outflows 
are between 0.3 x 10 5 yr and 6.1 x 10 5 yr. 85% 
of outflows have dynamical timescale shorter than 
2.5 x 10 5 yr. 

The mass, momentum, energy and luminosity 
in Table [5] are only lower limits because we did 
not take into account the inclination and blending 
correction in the calculation. The mass should be 
multiplied by a factor of 2 due to blending. As¬ 
suming the average inclination angle of outflows is 
57.3°, the velocity and the dynamic age should be 
scaled up by a factor of 1.9 and 0.64, respectively. 
Combining the correction factors due to blending 
and inclination, the momentum, the kinetic energy 
and luminosity of outflows should be multiplied by 
a factor of 3.8, 6.8 and 11, respectively. After 
correction, the total mass, momentum, energy, 
and luminosity of all outflows found in Taurus are 
approximately 15.4 Mq , 77 Mq km s -1 , 3.9 x 
10 45 erg, and 1.3 x 10 33 erg s 4 , respectively. 
The totals of the previously known outflows are 
about 8.6 Mq, 47 M 0 km s _1 , 2.6 x 10 45 
erg, and 9.4 x 10 32 erg s _1 , respectively. We 
found 1.8 times more outflowing mass, 1.6 times 
more momentum and 1.5 times more energy from 
outflows injecting into the Taurus molecular cloud 
than previously study. A high spatial dynamic 
range and systematic spectral line survey with 
good angular resolution is clearly necessary for 
obtaining a more complete picture of the influence 
of outflows on their parent cloud. 


4. Bubbles 


Following the method of identifying bubbles in 
Arce et al . ( 201 il l we have identified 37 bubbles in 
the ~ 100 deg 2 region of Taurus. The procedures 
for bubble searching, the morphology and physical 
parameters of bubbles are described in the follow¬ 


ing sections. 


4.1. The Procedures of Searching for Bub¬ 
bles 


We undertook a blind search for bubbles using 
the FCRAO 13 CO data cube. The integrated 
intensity map, P-V diagram and channel maps of 
each bubble were examined. The detailed steps of 
the search were as follows. 


(1) We first searched for circular or arc-like 
(hereafter bubble-like) structures in 13 CO data 
cube channel by channel through visual inspec¬ 
tion. If there is a bubble-like structure in at least 
three contiguous channels, we considered it as a 
bubble candidate. The approximate central posi¬ 
tion and radius of each candidate were recorded 
for further analysis. We also marked the channels 
where the bubble-like structure appears. With 
the marked channels we obtained the expanding 
velocity interval of a bubble. 


(2) We plotted 13 CO contour maps around the 
central position of the bubble candidates with the 
expanding velocity intervals. 

(3) We plotted P-V diagrams in 13 CO through 
the central position of each candidate of every 15° 
in position angle. We chose the one with the most 
obvious circular or “V” structure to show in the 
figures. The circular or “V” structure in the P- 
V diagram is desc ribed in the expanding bubble 


model (jArce et al. 2011). 


(4) We plotted the 13 CO channel maps of each 
candidate to look over the variation of radius with 
velocity. 

(5) Finally, we fitted a Gaussian profile to the 
azimuthally averaged profile of 13 CO intensity of 
each candidate in the channel where the bubble 
morphology is most like a ring or arc. The radius 
of a bubble was obtained from the peak position 
of the fitted profile. 


The contour map, P-V diagram, channel maps 
and Gaussian fitted curves helped us not only to 
analyze the morphology but also to determine the 
confidence level of a bubble. The bubble candi¬ 
dates were classified into 6 categories according to 
the characteristic of the above four types of plots. 
The criteria for bubble classification, as well as the 
numbers and ratios of bubbles in different classes 
are illustrated in Table [6] In this table “x” means 
that it meets a certain condition. For each type of 
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Table 5 —Continued 


Outflow 

Name 

Lobe 

^avg 

(km s _1 ) 

Area b 

(arcmin) 

Length 

(pc) 

Mass 

(M 0 ) 

Momentum 
(Mq km s ~ 1 ) 

Energy 
(10 43 erg) 

^•dyn 

(10 5 yr) 

-^flow 

(10 30 erg s -1 ) 

TMO-30 

Blueshifted 

_ 

_ 

_ 

_ 

_ 

_ 

_ 

_ 


Redshifted 

2.4 

9x9 

0.49 

0.286 

0.673 

1.57 

2.1 

2.43 

TMO.31 

Blueshifted 

- 

- 

- 

- 

- 

- 

- 

- 


Redshifted 

2.1 

5 X 14 

0.62 

0.277 

0.587 

1.24 

2.9 

1.37 

TMO.32 

Blueshifted 

- 

- 

- 

- 

- 

- 

- 

- 


Redshifted 

3.3 

9x6 

0.42 

0.131 

0.429 

1.39 

1.2 

3.53 

TMO.33 

Blueshifted 

2.1 

8x5 

0.38 

0.120 

0.249 

0.51 

1.8 

0.90 


Redshifted 

3.7 

11 x 13 

0.68 

0.896 

3.299 

12.08 

1.8 

1.174 

TMO.34 

Blueshifted 

3.2 

5x3 

0.23 

0.061 

0.198 

0.63 

0.7 

2.86 


Redshifted 

- 

- 

- 

- 

- 

- 

- 

- 

TMO_35 

Blueshifted 

- 

- 

- 

- 

- 

- 

- 

- 


Redshifted 

3.2 

6x4 

0.28 

0.049 

0.156 

0.49 

0.9 

1.79 

TMO.36 

Blueshifted 

2.6 

9 X 14 

0.65 

0.484 

1.241 

3.17 

2.5 

4.03 


Redshifted 

1.9 

8x6 

0.42 

0.092 

0.177 

0.34 

2.1 

0.51 

TMO.37 

Blueshifted 

- 

- 

- 

- 

- 

- 

- 

- 


Redshifted 

2.1 

15 x 19 

0.99 

0.348 

0.737 

1.55 

4.6 

1.08 

TMO.38 

Blueshifted 

2.3 

13 x 8 

0.61 

0.377 

0.860 

1.95 

2.6 

2.35 


Redshifted 

- 

- 

- 

- 

- 

- 

- 

- 

TMO.39 

Blueshifted 

2.3 

7x5 

0.35 

0.136 

0.313 

0.72 

1.5 

1.51 


Redshifted 

2.2 

6x5 

0.32 

0.076 

0.167 

0.37 

1.4 

0.81 

TMO.40 

Blueshifted 

- 

- 

- 

- 

- 

- 

- 

- 


Redshifted 

4.7 

4x4 

0.25 

0.129 

0.607 

2.84 

0.5 

7.683 

TMO.41 

Blueshifted 

3.8 

3x4 

0.23 

0.028 

0.104 

0.39 

0.6 

2.09 


Redshifted 

3.9 

5 x 15 

0.63 

0.317 

1.226 

4.72 

1.6 

9.401 

TMO.42 

Blueshifted 

3.2 

5 x 15 

0.65 

0.105 

0.338 

1.08 

2.0 

1.72 


Redshifted 

2.3 

3x3 

0.17 

0.004 

0.010 

0.02 

0.7 

0.10 

TMO.43 

Blueshifted 

2.7 

9x9 

0.51 

0.034 

0.092 

0.24 

1.9 

0.41 


Redshifted 

- 

- 

- 

- 

- 

- 

- 

- 

TMO.44 

Blueshifted 

2.4 

15 x 14 

0.84 

0.148 

0.359 

0.87 

3.4 

0.81 


Redshifted 

1.4 

9 x 13 

0.63 

0.070 

0.097 

0.13 

4.4 

0.10 

TMO.45 

Blueshifted 

2.9 

2x4 

0.19 

0.013 

0.037 

0.11 

0.6 

0.53 


Redshifted 

2.1 

3x2 

0.14 

0.006 

0.013 

0.03 

0.6 

0.14 

TMO.46 

Blueshifted 

2.8 

6 x 13 

0.58 

0.050 

0.141 

0.39 

2.0 

0.61 


Redshifted 

2.9 

5 x 18 

0.76 

0.123 

0.362 

1.06 

2.5 

1.34 

TMO.47 

Blueshifted 

1.5 

8x3 

0.36 

0.193 

0.282 

0.41 

2.4 

0.53 


Redshifted 

1.9 

5 x 10 

0.45 

0.304 

0.583 

1.11 

2.3 

1.53 

TMO-48 

Blueshifted 

5.2 

2x3 

0.15 

0.002 

0.012 

0.06 

0.3 

0.71 


Redshifted 

3.4 

2x2 

0.11 

0.001 

0.004 

0.02 

0.3 

0.16 

TMO.49 

Blueshifted 

- 

- 

- 

- 

- 

- 

- 

- 


Redshifted 

4.4 

7 X 17 

0.76 

0.042 

0.182 

0.79 

1.7 

1.47 

TMO.50 

Blueshifted 

- 

- 

- 

- 

- 

- 

- 

- 


Redshifted 

2.2 

4 x 13 

0.55 

0.020 

0.044 

0.10 

2.4 

0.12 

TMO.51 

Blueshifted 

- 

- 

- 

- 

- 

- 

- 

- 


Redshifted 

2.9 

2x2 

0.10 

0.001 

0.003 

0.01 

0.4 

0.08 

TMO.52 

Blueshifted 

- 

- 

- 

- 

- 

- 

- 

- 


Redshifted 

2.3 

2x3 

0.16 

0.002 

0.004 

0.01 

0.7 

0.04 

TMO.53 

Blueshifted 

- 

- 

- 

- 

- 

- 

- 

- 


Redshifted 

2.4 

2 X 10 

0.41 

0.005 

0.011 

0.03 

1.7 

0.05 

TMO.54 

Blueshifted 

1.9 

5x5 

0.30 

0.013 

0.025 

0.05 

1.5 

0.10 


Redshifted 

1.8 

8x5 

0.38 

0.034 

0.061 

0.11 

2.1 

0.17 

TMO.55 

Blueshifted 

- 

- 

- 

- 

- 

- 

- 

- 


Redshifted 

2.0 

4 x 11 

0.49 

0.011 

0.023 

0.05 

2.4 

0.06 


a The average velocity of the outflow relative to the cloud systemic velocity. 
b The extent along right ascension X the extent along declination. 
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Fig. 2.— Histograms of the distributions of outflow parameters. The shaded histograms represent the blue 
lobes and the open histograms represent the red lobes. 





plot the condition is as follows. 

(a) There is an obvious bubble-like structure in 
the contour map. 

(b) The P-V diagram has an obvious circular or 
“V” structure. 

(c) There is an obvious bubble-likc structure 
in the channel map and the radius of bubble is 
increasing or decreasing with channel. 

(d) The average intensity distribution can be 
fitted with a Gaussian profile. 

Meeting all the above four items is required for 
a high ranking (Class A). If the plots only meet (b) 
and (c), we then divided the bubbles into the lower 
ranking (Class Bl) because only the expanding 
velocity is detected but there is no obvious bubble¬ 
like structure and good Gaussian fitted profile. 
Then B2, B3 and B4 are in descending order of 
ranking. A candidate bubble only meeting (a) is 
assigned the lowest ranking (Class C) because the 
gas with expanding velocity is not obvious. 

Following the above procedures we found 37 
bubbles in the entire 100 deg 2 area of Taurus 
molecular cloud. Each bubble is referred as a 
“Taurus Molecular Bubble” (TMB). The positions 
and classifications are listed in Table 0 The num¬ 
bers and percentage of each class are illustrated in 
Table H 

4.2. Morphology of Bubbles 

The 13 CO integrated intensity maps, P-V dia¬ 
grams, Gaussian fitting profiles and channel maps 
for the bubbles are presented in Fig. [M] - Fig. 
[Ml If the morphology of a contour map is 
a closed ring, we then considered it to be an 
expanding bubble. If the ring on the contour map 
is incomplete, we then called it a broken bubble. 
There are 3 expanding bubbles (TMB_07, TMB_10 
and TMB_24) and 34 broken bubbles among all 
the bubbles in Taurus. 

4.3. Physical Parameters of Bubbles 

To examine the impact of bubbles on the host 
cloud we calculated the mass, momentum, kinetic 
energy, dynamical timescale and energy deposition 
rate of the bubbles. Assuming the 13 CO(1-0) 
emission of the bubble is optically thin, the total 


column density is derived as follows: 


3 k 2 T ex f T 


N ^ C °) = ^^eM-h»/kT ex) J T -*>' 

( 3 ) 

where k = 1.38 x 10 16 erg/K, h = 6.626 x 10 27 
erg s, /Ltd = 0.112 x 10~ 18 esu and u = 110.2014 x 
10 9 Hz. The excitation temperature, T ex is 
assumed to be 25 K. T s is the observed source 
antenna temperature with proper correction for 
antenna efficiency. The optical depth correction 
factor, f T , is esti mated from the following formulae 
( Qian et aI1l2012t) . 


fr = 


J r( 13 CO)dv 
/[I — e~' r ( 13 CO)]ch; ’ 


( 4 ) 


where r( 13 CO) is the opacity of the 13 CO transi¬ 
tion. Assuming equal excitation temperatures for 
13 CO and 12 CO, we can get 


t( 12 co) i - e - r ( 12c °) 

T( 13 CO) ~ 1 - e - T ( 13c °) ’ 


( 5 ) 


where T( 12 CO) and T( 13 CO) are the bright¬ 
ness temperature of 12 CO and 13 CO, respectively. 
r( 12 CO) is the opacity of the 12 CO transition. 
Assuming the 12 CO emission from the bubbles 
is optically thick (t( 12 CO) 1), the opacity of 

13 CO can be obtained from 


r( i3 CO) = — In I 1 — 


T( 13 CO) 

T( 12 CO) 


( 6 ) 


With the column density and area we can obtain 
the bubble mass. Using the bubble mass and 
expansion velocity we can then get the momentum 
and kinetic energy of the bubble using Pbubbie = 
Afbubbichexp and /’-.'bubble 2 -^bubble^exp> respec¬ 

tively. The kinetic timescale of bubble can be 
calculated as /kinetic = -R/Kxpj where R is the 
radius and I4xp is the expansion velocity of the 
bubble. The bubble energy injection rate, /.bubble; 
can be estimated as /.bubble — /’.'bubb ie// kinetic- 
The physical parameters of all bubbles are 
listed in Table 0 The momentum and kinetic en¬ 
ergy are lower limits mainly because of the under¬ 
estimate of the minimum expansion velocity. The 
total mass, momentum, energy and energy injec¬ 
tion rate of all detected bubbles in Taurus molec¬ 
ular cloud are about 1704 M 0 , 3780 M 0 km s _1 , 
9.2 x 10 46 erg and 6.4 x 10 33 erg s _1 , respectively. 


15 










Table 6 

Criteria for Bubble Classification and Bubble Distribution 


Bubble 

Class 

Contour 

Map 

P-V 

Diagram 

Channel 

Maps 

Fitting 

Curve 

Bubble 

Numbers 

Percentage 

A 

X 

X 

X 

X 

13 

35.2% 

B1 


X 

X 


6 

16.2% 

B2 

X 


X 


4 

10.8% 

B3 

X 

X 



4 

10.8% 

B4 



X 


4 

10.8% 

c 

X 




6 

16.2% 
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The distribution of radius, mass, energy and 
dynamical timescale of bubbles are shown in Fig. 
[3l The radius of bubbles is in the range 0.28 - 1.9 
pc. 78% of the bubbles are smaller than 1 pc. The 
mass of 65% of the bubbles is between 5 M 0 and 
50 M 0 . The bubbles with mass lower than 5 Mg 
and higher than 50 M 0 account for 16% and 19%, 
respectively. The highest bubble mass is 386 Mq. 
The energy of 60% of bubbles is in the range 
10 44 - 10 45 erg. The bubbles with energy lower 
than 10 44 erg and higher than 10 45 erg account 
for 16% and 24%, respectively. The dynamical 
timescales of bubbles are between 10 5 yr and 10 6 
yr. Almost 95% of bubbles are younger than 10 6 
yr. Compared to the outflows, the bubbles have 
about 110 times larger mass and 24 times higher 
energy. The extents of bubbles are larger than 
outflows, which can be seen from Fig. [4] The 
dynamical timescales of bubbles are longer than 
that of outflows. 


5. Analysis and Discussion 


5.1. Polarity of Outflows in Taurus 

Among the 55 outflows we found that bipolar, 
monopolar redshifted and monopolar blueshifted 
outflows account for 45%, 44% and 11%, respec¬ 
tively. There are more red lobes than blue ones , 
which can be seen from the histograms in Fig. [2] 
The occurrence of more red lob es may res ult from 


the fact that Taurus is thin ( Qian et ah 2014). 


Red lobes tend to be smaller and younger. The 
total mass and energy of red lobes is similar to 
blue lobes on average, which can be seen from the 
upper right panel and lower left panel of Fig. [2] 


5.2. The Driving Sources of Outflows and 
Bubbles in Taurus 

The outflows are driven by four types of YSOs. 
From Table [I] we can see Class I, Flat, Class II 
and Class III account for 36.4%, 18.2%, 27.3% 
and 12.7% of all the driving sources, respectively. 
Fig. [Ushows the distribution of different classes of 
YSOs driving outflows (hereafter outflow-driving 
YSO) and YSOs inside the bubbles (hereafter 
bubble-driving YSO). The rough dividing line 
shows that there are more outflow-driving YSOs 
in Class I, Flat and Class II while few outflow¬ 
driving YSOs in Class III, which indicates that 
outflows are more likely appear in the earlier stage 
(Class I) than in the later phase (Class III) of 
star formation. There are more bubble-driving 
YSOs of Class II and Class III while there are few 
bubble-driving YSOs of Class I and Flat, implying 
that the bubble structures are more likely to occur 
in the later stage of star formation. From the size 
of the symbols we can see that the larger outflows 
and bubbles are, the higher energy they have. 

5.3. The Feedback of Outflows and Bub¬ 
bles in Taurus 

With the complete sample of outflows and 
bubbles we can estimate the overall impact of 
dynamical structures on the Taurus molecular 
cloud. We investigated whether the outflows and 
bubbles have enough energy to potentially unbind 
the entire Taurus molecular cloud or drive the 
turbulence in the cloud. 


5 . 3 . 1 . The Energy of Outflows and Bubbles 
Cannot Balance the Gravitational 
Binding Energy of the Entire Cloud 


Using a total mass of 1.5 xlO 4 M 0 (Pineda^t_al 


2010) and an effective radius of 13.8 pc (INaravanan et al 
20121) for the 100 deg 2 region of Taurus, we 


calculated the magnitude of the gravitational 
binding energy GM c 2 loud /i? c i oud to be ~ 1.5 x 10 48 
erg. The total kinetic energy of outflows from the 
44 deg 2 region of Taurus is 3.9 x 10 45 erg, much 
less than the gravitational binding energy. Given 
that we searched for outflows around YSOs only in 
the Spitzer 44 deg 2 survey region not the overall 
area of Taurus, we may well have missed some 
outflows. Most of the gas of Taurus is centered on 
the Spitzer 44 deg 2 survey region, which can be 


17 













Table 7 

BUBBLES IN TAURUS 


Bubble 

Name 

R.A. 

(J2000) 

Dec. 

(J2000) 

Bubble 

Class 

YSO 

Radius 

(pc) 

V exp 

(km s -1 ) 

Mass 

(M 0 ) 

Momentum 
(Mq km s -1 ) 

Energy 
(10 45 erg) 

^■kinetic 

10 6 yr 

^bubble 

(10 32 erg s -1 ) 

m w a 

(10- 8 Mq yr" 1 ) 

TMB.Ol 

04 12 08 

24 53 33 

A 

N 

0.98 

1.3 

25 

31 

0.39 

0.8 

0.17 

15.5 

TMB.02 

04 14 28 

27 45 53 

B1 

Y 

0.60 

1.3 

6 

7 

0.09 

0.5 

0.06 

3.5 

TMB.03 

04 16 20 

28 28 53 

A 

Y 

0.76 

1.8 

10 

18 

0.33 

0.4 

0.25 

9.0 

TMB.04 

04 19 05 

27 33 33 

A 

Y 

0.62 

1.3 

7 

9 

0.12 

0.5 

0.08 

4.5 

TMB.05 

04 21 12 

26 55 33 

B3 

Y 

0.56 

1.8 

4 

8 

0.14 

0.3 

0.14 

4.0 

TMB.06 

04 25 17 

25 32 13 

B1 

N 

0.77 

1.3 

10 

12 

0.16 

0.6 

0.08 

6.0 

TMB.07 

04 25 29 

26 10 13 

B1 

N 

1.12 

2.3 

102 

234 

5.31 

0.5 

3.51 

117.0 

TMB.08 

04 27 07 

24 20 13 

A 

N 

0.70 

2.3 

18 

42 

0.95 

0.3 

1.00 

21.0 

TMB.09 

04 27 31 

26 16 53 

B4 

N 

0.49 

1.0 

4 

4 

0.04 

0.5 

0.03 

2.0 

TMB.10 

04 28 52 

24 14 33 

A 

Y 

1.58 

1.5 

213 

325 

4.92 

1.0 

1.53 

162.5 

TMB.ll 

04 29 32 

26 32 33 

A 

Y 

0.84 

2.0 

41 

83 

1.68 

0.4 

1.31 

41.5 

TMB.12 

04 29 44 

26 32 53 

B4 

Y 

0.70 

2.5 

23 

58 

1.47 

0.3 

1.72 

29.0 

TMB.13 

04 30 31 

24 26 13 

A 

Y 

0.73 

2.3 

18 

41 

0.93 

0.3 

0.94 

20.5 

TMB.14 

04 31 14 

29 25 53 

B2 

Y 

0.84 

2.0 

10 

20 

0.40 

0.4 

0.31 

10.0 

TMB.15 

04 31 30 

24 14 33 

A 

Y 

0.84 

1.3 

34 

43 

0.54 

0.6 

0.27 

21.5 

TMB.16 

04 31 32 

24 09 53 

B3 

N 

0.70 

1.8 

22 

39 

0.69 

0.4 

0.57 

19.5 

TMB.17 

04 31 35 

23 35 13 

B4 

N 

0.42 

1.0 

2 

2 

0.02 

0.4 

0.02 

1.0 

TMB.18 

04 31 50 

24 22 13 

c 

N 

0.56 

2.0 

14 

29 

0.59 

0.3 

0.69 

14.5 

TMB.19 

04 31 59 

25 43 13 

B1 

Y 

0.70 

2.0 

16 

32 

0.64 

0.3 

0.60 

16.0 

TMB.20 

04 32 03 

25 36 53 

B2 

N 

0.42 

1.3 

4 

5 

0.07 

0.3 

0.07 

2.5 

TMB.21 

04 32 37 

29 29 13 

A 

N 

0.62 

2.3 

6 

14 

0.31 

0.3 

0.37 

7.0 

TMB.22 

04 32 39 

24 46 13 

A 

Y 

1.06 

1.3 

28 

36 

0.45 

0.8 

0.17 

18.0 

TMB.23 

04 33 10 

26 08 53 

B4 

N 

0.28 

1.3 

2 

2 

0.03 

0.2 

0.04 

1.0 

TMB.24 

04 33 13 

25 24 53 

c 

N 

0.49 

1.3 

5 

7 

0.09 

0.4 

0.07 

3.5 

TMB.25 

04 33 34 

24 20 53 

A 

Y 

0.28 

3.0 

5 

15 

0.46 

0.1 

1.61 

7.5 

TMB.26 

04 34 47 

29 37 13 

B3 

N 

0.46 

2.3 

3 

6 

0.14 

0.2 

0.23 

3.0 

TMB.27 

04 36 02 

28 23 13 

C 

Y 

1.40 

2.5 

64 

161 

4.08 

0.5 

2.39 

80.5 

TMB.28 

04 36 23 

25 36 33 

B2 

Y 

1.40 

3.3 

386 

1275 

41.84 

0.4 

31.93 

637.5 

TMB.29 

04 37 04 

25 46 33 

c 

Y 

0.84 

1.5 

20 

30 

0.46 

0.5 

0.27 

15.0 

TMB.30 

04 38 11 

26 05 53 

B3 

Y 

1.90 

1.8 

340 

604 

10.68 

1.0 

3.23 

302.0 

TMB.31 

04 39 11 

29 05 13 

B2 

N 

1.26 

1.5 

74 

113 

1.71 

0.8 

0.67 

56.5 

TMB.32 

04 39 48 

28 35 33 

c 

N 

0.70 

2.0 

7 

14 

0.27 

0.3 

0.26 

7.0 

TMB.33 

04 41 10 

25 31 13 

c 

N 

0.70 

1.0 

10 

10 

0.10 

0.7 

0.05 

5.0 

TMB.34 

04 44 20 

28 36 53 

B1 

N 

1.26 

2.8 

143 

399 

11.08 

0.4 

7.95 

199.5 

TMB.35 

04 46 12 

25 07 33 

A 

N 

0.62 

1.5 

8 

13 

0.19 

0.4 

0.15 

6.5 

TMB.36 

04 46 43 

24 59 13 

B1 

Y 

0.56 

1.8 

12 

21 

0.38 

0.3 

0.39 

10.5 

TMB.37 

04 48 12 

24 50 33 

A 

N 

0.63 

2.3 

8 

18 

0.41 

0.3 

0.48 

9.0 


'Estimate of minimum stellar wind mass loss rate needed to drive the bubbles. 
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Fig. 4.- Distribution of Outflows and Bubbles. The open upward triangles represent the “outflow-driving 
YSOs” and the open circles represent the “bubble-driving YSOs”. The black solid line roughly divides the 
outflows and bubbles. The size of symbols is proportional to the energy of outflows and bubbles. 








seen from Fig. 2 of Goldsmith et al.l ( 2008h . There 
are not many YSOs outside the Spitzer coverage 
in Taurus and those YSOs are generally clustered, 
which can be seen from Figure 1 of lRebull et al. 
(201111. So there should be few outflows outside 
the Spitzer coverage in Taurus and the outflows we 
found around YSOs in the 44 deg 2 area account for 
the majority of outflows in Taurus. Similarly, the 
total kinetic energy of the detected bubbles in the 
100 deg 2 Taurus region is 9.2 x 10 46 erg, which also 
cannot balance the gravitational potential energy 
of the entire cloud. 


5.3.2. Turbulent Energy is Greater than 
the Energy of Outflows and Bubbles 

The turbulent energy of the Taurus 
molecular cloud is given approximately by 

Eturb = -McloudOM, (7) 


where 03 d is the three dimensional turbulent 
velocity dispersion, which can be calculated 
by 

03d = — Av FWH m- (8) 

2\/2 In 2 

Here Avfwhm = 2 km s _1 is the one 
dimensional FWHM velocity dispersion 
based on typical 13 CO spectra in Taurus 
( Narayanan et al. 2012 ). Then we get 173 d = 
1.47 km s _1 . The total mass of the 100 deg 2 
of Taurus is M c i ou d = 1.5 x 10 4 M,r 

2010h . 


region 


(Pineda et al. 


J-cloud 

Using Eq. 


m 


© 

we 

obtain the turbulent energy of the Taurus 

to be 3.2 x 10 47 erg. The energy of all detected 
outflows (3.9 x 10 45 erg) is about two orders of 
magnitude less than the turbulent energy of the 
cloud. The lower limit of the total energy of the 
bubbles (9.2 x 10 46 erg) is 29% of the turbulent 
energy. We conclude that the total energy 
of outflows and bubbles cannot balance the 
turbulence in Taurus. 


5.3.3. Turbulent Dissipation Rate is Com¬ 
parable to Bubble Energy Injection 
Rate 

We also estimated the total outflow energy rate 
(outflow luminosity, Tfl ow ) and the total bubble 
energy rate (bubble luminosity, Lbubbie) with the 
energy rate needed to maintain the turbulence 


(turbulent energy dissipation, L tU rb)- The lu¬ 
minosity of outflows is 1.3 x 10 33 erg s _1 
after inclination and blending correction. 
Summing up the luminosity in Table [5] we 
get the energy injection rate of bubbles to 
be 6.4 x 10 33 erg s” 1 . 

The turbulent dissipation rate can be 
calculated as 


Eturb — 


E- 


turb 


f dif 


(9) 


where tdi ss is the turbulent dissipation time. 
We estimate the turbulent dissipation time 
through two methods based on numerical 
simulations. 

First, the turbulent dissipation time of 
the cloud is given by (McKee &: Ostriker 
20071 ) 

t diss ~ 0.5 — , (10) 

0’ld 

where d = 27.6 pc is the cloud diameter 
and (7 id is the one dimensional turbulent 
velocity dispersion along the line of sight, 


Av 


0’ld = 


FWHM 


2 V 2 In 2 


( 11 ) 


Here Avfwhm = 2 km s _1 is the same as 
that in § 15. 3.21 Then we get erid = 0.85 
km s _1 . Combining Eq. (11011 and Eq. 
m we obtain the turbulent dissipation 
time, tdiss = 1.6 x 10 7 yr, which is about 
6 times_Jarger than the result (2.7 x 10 6 yr) 
in iNaravanan et al. ( 2012! ). Then using Eq. 
(l9l) we get the turbulent dissipation rate to 
be 6.6 x 10 32 erg s -1 , which is 51% of the 
luminosity of outflows and only 10 % of the 
energy injection rate of bubbles. 

Second, we follow Mac Low ( 19991 ) to 
calculate the dissipation time of the cloud 
by 

/ 3.9k 

tdiss ~ (tt )tff, (12) 

-LVJ-rms 

where tff is the free-fall timescale, M rms is 
the Mach number of the turbulence, and 
k is the ratio of the driving length to the 
Jean’s length of the cloud, 


Ad 


(13) 
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Using our extensive data sets we get Ad = 
0.53 pc, which is the average size of outflows 
and bubbles we found in Taurus. The Jean’s 
length is defined as 


Aj = c s (7r/Gp reg ) 1/2 , 


(14) 


wher e c s is the sound speed ( Mac Low! 

(15) 


1999). For an ideal gas 

c s = (3kT/^ 2 , 


where k is the Boltzmann’s constant, T 
= 10 K is the temperature of the Taurus 
molecular cloud, and /i = 2.72 is the mean 
molecular weight ( Brunt 2010). Then we 
get the sound speed, c s = 0.3 km s -1 . p Teg 
is the representative volume density of the 
region where dissipation takes place. We 
estimate the volume density to be ~ 1500 
cm -3 . Using Eq. (114(1 we obtain the Jean’s 
length of the region, Aj = 0.83 pc, which 
is about 4 time s larger than th at (Aj ~ 0.2 
pc) in Perseus (Arce et al. 20100 . And then 
using Eq. (113ft we obtain n = 0.64, which 
is differ ent fr o m th e assumpt ion (k ~ 1) 
by Arce et al. ( 2010p and Narayanan et al. 

(<2012h . 


The Mach numb er of the turb ulence can 
be calculated by ( Mac Low 19991) 


& 3d 


(16) 


Using Eq. ([8]), Eq. (1151) and Eq. (116ft we 
get the Mach number, M rms = 5, which is 
dif ferent from the a ssumption (M r rns = 10) 
by Arce et al. ( 2010p and Narayanan et al. 

(|2012h . 

The free-fall timescale of the cloud. 


t ff = (37T/32G/9, 


cloud 


, 1/2 


(17) 


where 


Pcloud — 


3M c ioud 


47rR cloud 


(18) 


is the average volume density of the cloud. 
Then we get the free-fall timescale, tff = 
7 x 10 6 yr. 

Using the formulas from Eq. (1121) to 
Eq. (1181) we obtain the turbulent dis¬ 
sipation time to be 3.5 x 10 6 yr. Then 


we get the turbulent dissipation rate to 
be 3.1 x 10 33 erg s -1 , which is about 2.4 
times larger than the luminosity of outflows 
but 48% of the energy injection rate of 
bubbles. The turbulent dissipation rate 
we obtained _is close to that (3.8 x 10 33 
erg s _1 ) of Narayanan et ah] (12012). In 
Table [8] we list the parameters related to 
the dissipation rate we get from the above 
two methods and compare them with the 
results of Narayanan et al. (2012). 

Both methods invoke numerical simula¬ 
tions to calibrate the numerical factors in 
addition to essentially dimensional argu¬ 
ments. The main difference of the two 
methods is the scale of the region where 
dissipation takes place. McKee & Ostriker 
(120071) adopted the dimension of the en¬ 
tire cloud, while when we use the method 


given by Mac Low ( 1999f ) the scale is the 
average size of outflows and bubbles. None 
of the simulations so far implements the 
physics (excitation, radiative transfer, etc.) 
necessary for actually modeling dissipation. 
Thus we should treat the calculations above 
with caution and take them as dimensional 
and order-of magnitude estimates. 

Comparing the energy injection rate of 
outflows and bubble with the turbulent 
dissipation rate, we conclude that in the 
current episode of star formation in Taurus, 
both outflows and bubbles can sustain the 
currently observed turbulence in Taurus. 


5.4. Protostellar Winds Can Drive Bub¬ 
bles to Sustain Turbulence in Taurus 


Protostellar winds 
the cloud and may 


will inject energy into 
help sustain turbulence 
([Nakamura & Li 20071 ). The winds can clear the 
gas surrou ndin g the young star and form a bubble 
structure ( Arce et aDl2011 ). To assess whether the 
winds can drive bubbles in Taurus we compared 
the wind energy injection rate into the cloud (E w ) 
with the total energy injection rate f rom bubbles 
of the cloud. Following lArce et ah ( 201ll ). we 
estimated the wind_ e nerg y injection rate using 
Equation 3.7 from [McKee (1989): 


E w — ( M W V W )u” 3 d 5 


(19) 
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where v w is the wind velocity, which is generally 
assumed to be close to the star escape velocity. 
For the low- and intermediate-mass stars the 
escape v elocit y is a bout 1 — 4 x 10 2 km s -1 


( Lamers fc Cassinelli 19991. Similar to Arce et al 
( 20111) we assumed v w ~ 200 km s -1 . The total 


mass loss rate from the protostellar winds is given 
by M w , which can be estimated by the sum of 
the wind mass loss rate for each bubble (m w ). 
The wind mass loss rate required to produce the 
bubb les is roughly estimated by Equation (2) from 
Arce et al. ( 2011 1: 


m w = 


(bubble 




( 20 ) 


where Pbubbie is the total momentum of bubbles. 
The wind velocity, v w , is the same as that in Eq. 
m- r w is the wind timescale, which is assumed 


1 Myr ( Arce et al. 120111) . From Eq. (l20l) we 
obtained the wind mass loss rates of each bubble 
(to w ), which are listed in Table |T| Summing 
m w of all bubbles we find M w to be 1.89 x 10 -5 
Mg yr -1 . Using Eq. (flTH) we obtained the wind 
energy injection rate (U w ) to be ~ 2 x 10 33 erg 
s -1 , 31% of the total energy injection rate from 
bubbles in Taurus, which is comparable to the 
turbulent dissipation rate in Taurus. Therefore, 
the protostellar winds can drive bubbles to sustain 
turbulence in Taurus. 


5.5. Potential Sources of Turbulent Mo¬ 
tions in Taurus 

The origin of turbulence in the molecular 
cloud has been intensely debated over 
the past three decades (e.g. 


Larson 


1981; 


Hev er &: Bruntll2004f ). iHennebelle &: Falgarone 
(|2012h suggests that for a large frac¬ 
tion of clouds the turbulent driving is 
external. Numerical simulations shows 
that the external sources of turbulence 
likely to be large-scale Hi streams 


are 


(Ballesteros-Paredes et al. 1999h . shocks 

"HL. 


(Nakamura & Li 

2007 


Wane et al. 2010), 

supernovae explosion 

and galactic dif- 

ferential rotation ( 

Klessen & Hennebelle 


2010; iHennebelle &; Falgarone 2012). 


is unclear which one 
turbulence in Taurus. 


is the source 


It 

of 


6. Conclusions 

We have studied the dynamic structures in¬ 
cluding outflows and bubbles within the Taurus 
molecular cloud using the 100 deg 2 FCRAO large- 
scale 12 CO(1-0) and 13 CO(1-0) maps and the 
Spitzer protostellar catalog. The high sensitivity 
and large spatial dynamic range of the maps 
provide us an excellent opportunity to undertake 
an unbiased search for outflows and bubbles in 
this region. We also analyzed the energy injection 
of these dynamic structures into the entire cloud. 
Our conclusions regarding the dynamic structures 
in Taurus and their properties are as follows. 

1. We identified 55 outflows around the Spitzer 
YSOs in the main 44 deg 2 area of Taurus. In 
total, 31 of the detected outflows were previously 
unknown, increasing the number of outflows by a 
factor of 1.3. 

2. We classified the outflows into 5 categories 
according to the morphology of contour maps and 
P-V diagrams. The classifications indicate the 
confidence level of the outflows. 76.3% of the 
outflows are in the “most probable” category in 
our study. 

3. Most of the outflows are driven by Class I, 
Flat and Class II YSOs while few outflows were 
found around Class III YSOs, which indicates 
that the outflow activity likely occurred in the 
earlier stage rather than the late phase of the star 
formation. 

4. More bipolar and monopolar redshifted 
outflows were identified while few monopolar 
blueshifted ones were detected in our study. 

5. We detected 37 bubbles in the 100 deg 2 
region of Taurus. All the bubbles were previously 
unknown. The bubbles were identified by the in¬ 
tegrated intensity maps, P-V diagrams, Gaussian 
fitting profiles and channel maps. 

6. The gravitational binding energy of the 
Taurus molecular cloud is ~ 1.5 x 10 48 erg. The 
total kinetic energy of outflows and bubbles in 
Taurus are ~ 3.9 x 10 45 erg and ~ 9.2 x 10 46 
erg, respectively. Neither outflows nor bubbles can 
balance the overall gravitational binding energy of 
Taurus. 

7. The turbulent energy of the Taurus molec¬ 
ular cloud is ~ 3.2 x 10 47 erg. The energy of 
all detected outflows and bubbles cannot have 
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generated the observed turbulence in Taurus. 

8. The rate of turbulent dissipation in Taurus is 

between ~ 6.6 x 10 32 to ~ 3.1 x 10 33 erg s _1 . 
The energy injection rates of outflows and bubbles 
are ~ 1.3 x 10 33 erg s _1 and ~ 6.4 x 10 33 erg 
s~ 1 . respectively. Both outflows and bubbles 
can sustain the turbulence in Taurus at the 
current epoch. 

9. The stellar winds can drive bubbles to 
sustain turbulence in the Taurus molecular cloud. 


Table 8: The Parameters Related to the Dissipa¬ 
tion Rate 


Method 

K 

^rms 

^diss 

(10 6 yr) 

Tturb 

(10 33 erg s 1 ) 

MO07 a 

- 

- 

16 

0.7 

ML99 b 

0.64 

5 

3.5 

3.1 

Nal2 c 

1 

10 

2.7 

3.8 


The resu l t we o bt ain u sing the method given by 


McKee & Ostriker ( 2007 ). b The result we obtain 


using th e method given by Mac Low ( 1999h . c The 
result oflNara^anaiiet_aL ( 20121) using the method 
given by Mac Low ( 19991) . 
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Fig. 5.— TMO_02 (SST 041412.2+280837): upper left, 12 CO integrated intensity map overlaid on the 
13 CO grey-scale image, integrated over 2 to 5 km s _1 for the blue lobe and 8 to 12.5 km s -1 for the red 
lobe with 12 CO, and integrated over 5.5 to 7.5 km s -1 for the 13 CO grey-scale image. The blue and red 
contour levels are 40%, 50%,..., 90% and 30%, 45%,..., 90% of their peak value, respectively. In this panel 
and all subsequent upper left panels including those of bubbles, the green filled stars, if present, show the 
location of Class I YSOs in Taurus, the green open squares show the location of flat-spectrum YSOs, the 
green open triangles show the location of Class II YSOs and the green fille d squares sho w the location of 


Class III YSOs. All the above YSOs were listed in Table 6 and Table 7 of ( Rebull et al. 2010 ). The black 


solid line represents a cut for the P-V diagram shown in the upper right panel, and the amaranth filled circle 
indicates the origin on the Y-axis of the upper right panel. The filled circle in the lower left corner shows 
the beam, upper right, P-V diagram of 12 CO, through the slice shown in the upper left panel at a position 
angle of 25°. Contour levels are 0.4 to 1.6 K by 0.3 K, 2 to 5 K by 0.5 K. lower left, average spectra of 
12 CO emission (blue lines) and 13 CO emission (red lines) towards the blueshifted lobe shown in the upper 
left panel, lower right, average spectra of 12 CO emission (blue lines) and 13 CO emission (red lines) towards 
the redshifted lobe shown in the upper left panel. 


25 









































Fig. 6.- TMO-Ol (SST 041159.7+294236): upper left and upper right, the same as upper two panels of 
Fig. [5] except for the following. In the upper left panel the integrated intervals for red lobe and grey-scale 
image are 8.75 to 10 km s^ 1 and 6.5 to 7.5 km s _1 , respectively. In the upper right panel the position angle 
is 20° and contour levels are 0.4 to 1.6 K by 0.3 K, 2 to 4 K by 0.5 K. lower, the same as the lower right 
panel of Fig. O 
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Fig. 7. TMO_03 (SST 041414.5+282758): The same as Fig. [5] except for the following. In the upper left 
panel the integrated intervals for red lobe and grey-scale image are 8.5 to 9.5 km s -1 and 6 to 7 km s _1 , 
respectively. The blue and red contour levels are 30%, 40%,..., 90% and 50%, 60%,..., 90% of their peak 
value, respectively. In the upper right panel the position angle is 45° and contour levels are 0.4 to 1.6 K by 
0.3 K, 2 to 6 K by 0.5 K. 


27 


Fig. 8 .- TMO_04 (SST 041832.0+283115): The same as Fig. [5] except for the following. In the upper left 
panel the integrated intervals for blue lobe, red lobe and grey-scale image are 4.4 to 5.3 km s -1 , 9.5 to 10.3 
km s _1 and 6.5 to 8 km s _1 , respectively. The red contour levels are 35%, 50%,..., %90 of the peak value. 
In the upper right panel the position angle is 0° and contour levels are 0.5 to 2 K by 0.3 K, 2.5 to 5.5 K by 
0.5 K. 
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Fig. 9.- TMOJ35 (SST 041941.4+271607): upper left and upper right, the same as upper two panels 
of Fig. [5] except for the following. In the upper left panel the integrated intervals for the blue lobe and 
grey-scale image are 2 to 4.3 km s -1 and 6 to 7 km s _1 , respectively. The blue levels are 30%, 45%,..., 90% 
of the peak value. In the upper right panel the position angle is 40° and contour levels are 0.3 to 0.9 K by 
0.15 K, 1.1 to 1.9 K by 0.2 K. lower, the same as the lower left panel of Fig. [5] 
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Fig. 10. TMCL06 (SST 041958.4+270957): The same as Fig. [5] except for the following. In the upper left 
panel the integrated intervals for the blue lobe, red lobe and grey-scale image are 2 to 4 km s _1 , 9 to 12 km 
s -1 and 5 to 7 km s _1 , respectively. The blue and red contour levels are both 30%, 40%,..., 90% of their 
peak value. In the upper right panel the position angle is 152° and contour levels are 0.5 to 2 K by 0.3 K, 
2.5 to 5.5 K by 0.5 K. 


30 


Fig. 11.— TMO_07 (SST 042107.9+270220): The same as Fig. [5] except for the following. In the upper 
left panel the integrated intervals for the blue lobe and red lobe are 3 to 4.5 km s -1 and 8 to 9 km s _1 , 
respectively. The blue and red contour levels are 50%, 60%,..., 90% and 40%, 50%,..., 90% of their peak 
value, respectively. In the upper right panel the position angle is 130° and contour levels are 0.5 to 2 K by 
0.3 K, 2.5 to 4.5 K by 0.5 K. 
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Fig. 12.— TMCL08 (SST 042215.6+265706): The same as Fig. [5] except for the following. In the upper left 
panel the integrated intervals for the blue lobe,red lobe and grey-scale image are 2 to 4.5 km s _1 , 8.3 to 9.5 
km s -1 and 5.5 to 7 km s _1 , respectively. The blue and red contour levels are both 30%, 40%,..., 90% of 
their peak value. In the upper right panel the position angle is 122° and contour levels are 0.5 to 2 K by 0.3 
K, 2.5 to 4 K by 0.5 K. 
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Fig. 13.- TMCL09 (SST 042325.9+250354): The same as Fig. [6] except for the following. In the upper left 
panel the integrated interval for the red lobe is 8.25 to 9.7 km s -1 . The red contour levels are 50%, 60%,..., 
90% of the peak value. In the upper right panel the position angle is 52° and contour levels are 0.4 to 1.6 K 
by 0.3 K, 2 to 4.5 by 0.5 K. 
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Fig. 14.— TMO_10 (SST 042420.9+263051): The same as Fig. [5] except for the following. In the upper left 
panel the integrated intervals for the blue lobe,red lobe and grey-scale image are 4 to 4.75 km s _1 , 7.5 to 
10 km s -1 and 5.5 to 7.25 km s _1 , respectively. The blue and red contour levels are 50%, 60%,..., 90% and 
60%, 70%,..., 90% of their peak value, respectively. In the upper right panel the position angle is 30° and 
contour levels are 0.4 to 1.6 K by 0.3 K, 2 to 7 K by 0.5 K. 
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Fig. 15.— TMCL11 (SST 042445.0+270144): The same as Fig. [6]cxcept for the following. In the upper left 
panel the integrated intervals for the red lobe and grey-scale image are 8 to 10.5 km s” 1 and 6 to 7 km s _1 , 
respectively. The red contour levels are 20%, 30%,..., 90% of the peak value. In the upper right panel the 
position angle is 0°. 
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Fig. 16.— TMCL12 (SST 042930.0+243955): The same as Fig. [6] except for the following. In the upper left 
panel the integrated intervals for the red lobe and grey-scale image are 8 to 9.2 km s -1 and 5.5 to 7.5 km 
s _1 , respectively. The red contour levels are 50%, 60%,..., 90% of the peak value. In the upper right panel 
the position angle is 53° and Contour levels are 0.4 to 4.8 K by 0.3 K. 
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Fig. 17.- TM0_13: upper left and upper right, the same as upper two panels of Fig. [6] except for the 

following. In the upper left panel the integrated intervals for the red lobe and grey-scale image are 7 to 9.5 
km s -1 and 5.5 to 6.5 km s _1 , respectively. The red contour levels are 40%, 50%,..., 90% of the peak value. 
In the upper right panel the position angle is 55° and contour levels are 0.4 to 1.6 K by 0.3 K, 2 to 4.5 K 
by 0.5 K. lower left, the same as the lower panel of Fig. [6] lower right, 13 CO integrated intensity map 
overlaid on the 13 CO grey-scale image, the integrated intervals for the red lobe and grey-scale image and the 
red contour levels are the same with those in the upper left panel of this figure. 
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Fig. 18.— TM0_14 (SST 043158.4+254329): The same as Fig. [5] except for the following. In the upper left 
panel the integrated intervals for the blue lobe,red lobe and grey-scale image are 3 to 4.5 km s _1 , 7 to 10 
km s -1 and 5 to 6.5 km s _1 , respectively. The blue and red contour levels are both 50%, 60%,..., 90% of 
their peak value. In the upper right panel the position angle is 145° and contour levels are 0.4 to 1.6 K by 
0.3 K, 2 to 4.5 K by 0.5 K. 
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Fig. 19.— TM0_15 (SST 043214.6+223742): The same as Fig. [5] except for the following. In the upper left 
panel the integrated intervals for the blue lobe,red lobe and grey-scale image are 2.8 to 4 km s _1 , 6.8 to 9 
km s -1 and 5 to 6 km s _1 , respectively. The blue and red contour levels are 60%, 70%,..., 90% and 40%, 
50%,..., 90% of their peak value, respectively. In the upper right panel the position angle is 55° and contour 
levels are 0.5 to 2 K by 0.3 K, 2.5 to 4 K by 0.5 K. 
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Fig. 20. TMO_16 (SST 043231.7+242002): The same as Fig. [5]except for the following. In the upper left 
panel the integrated intervals for the blue lobe,red lobe and grey-scale image are 3 to 4.5 km s _1 , 7.5 to 8.7 
km s -1 and 5 to 7 km s _1 , respectively. The blue and red contour levels are 80%, 90% and 50%, 60%,..., 
90% of their peak value, respectively. In the upper right panel the position angle is 20° and contour levels 
are 0.4 to 1.6 K by 0.3 K, 2 to 7 K by 0.5 K. 
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Fig. 21.— TMCL17 (SST 043232.0+225726): The same as Fig. [6]cxcept for the following. In the upper left 
panel the integrated intervals for the red lobe and grey-scale image are 6.8 to 7.8 km s -1 and 4.8 to 6.3 km 
s _1 , respectively. The red contour levels are 40%, 50%,..., 90% of the peak value. In the upper right panel 
the position angle is 170° and contour levels are 0.4 to 1.6 K by 0.3 K, 2 to 7 K by 0.5 K. 
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Fig. 22.— TMCL18 (SST 043243.0+255231): The same as Fig. [5] except for the following. In the upper left 
panel the integrated intervals for the blue lobe,red lobe and grey-scale image are 3.2 to 4 km s _1 , 7 to 8 
km s -1 and 5 to 6.5 km s _1 , respectively. The blue and red contour levels are 50%, 60%,..., 90% and 60%, 
70%,..., 90% of their peak value, respectively. In the upper right panel the position angle is 0° and contour 
levels are 0.4 to 1.6 K by 0.3 K, 2 to 4 K by 0.5 K. 
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Fig. 23.- TMCL19 (SST 043415.2+225030): The same as Fig. [6] except for the following. In the upper left 
panel the integrated intervals for the red lobe and grey-scale image are 7 to 8.5 km s -1 and 6 to 6.5 km s _1 , 
respectively. The red contour levels are 40%, 50%,..., 90% of the peak value. In the upper right panel the 
position angle is 125°. 
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Fig. 24.— TMO_20 (SST 043724.8+270919): The same as Fig. [6] except for the following. In the upper left 
panel the integrated intervals for the red lobe and grey-scale image are 8 to 10.5 km s -1 and 5 to 7 km s _1 , 
respectively. The red contour levels are 30%, 40%,..., 90% of the peak value. In the upper right panel the 
position angle is 90° and contour levels are 0.5 to 2 K by 0.3 K, 2.5 to 4 K by 0.5 K. 
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Fig. 25.— TM0_21 (SST 043953.9+260309): The same as Fig. [5] except for the following. In the upper left 
panel the integrated intervals for the blue lobe,red lobe and grey-scale image are 2 to 4 km s _1 , 8 to 10 km 
s -1 and 4.5 to 7.5 km s _1 , respectively. The blue and red contour levels are 30%, 45%,..., 90% and 40%, 
50%,..., 90% of their peak value, respectively. In the upper right panel the position angle is 90° and contour 
levels are 2 K by 0.3 K, 2.5 to 6 K by 0.5 K. 
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Fig. 26.- TMO_22 (SST 044108.2+255607): The same as Fig. [9] except for the following. In the upper left 
panel the integrated intervals for the blue lobe and grey-scale image are 3 to 4.5 km s -1 and 4.6 to 7 km 
s _1 , respectively. The blue contour levels are 35%, 50%,..., 95% of the peak value. In the upper right panel 
the contour levels are 0.4 to 1.6 K by 0.3 K, 2 to 4.5 K by 0.5 K. 
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Fig. 27. TMO_23 (SST 044112.6+254635): The same as Fig. [6]cxcept for the following. In the upper left 
panel the integrated intervals for the red lobe and grey-scale image are 7.5 to 10.5 km s -1 and 4.5 to 7 km 
s _1 , respectively. The red contour levels are 40%, 50%,..., 90% of the peak value. In the upper right panel 
the position angle is 0°. 
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Fig. 28.— TMO_24 (SST 044207.7+252311): The same as Fig. [6]cxcept for the following. In the upper left 
panel the integrated intervals for the red lobe and grey-scale image are 8 to 9.5 km s -1 and 5 to 7 km s _1 , 
respectively. The red contour levels are 30%, 40%,..., 90% of the peak value. In the upper right panel the 
position angle is 15° and contour levels are 0.4 to 1.6 K by 0.3 K, 2 to 4.5 K by 0.5 K. 
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Fig. 29. TMO_25 (SST 041858.1+281223): upper left and upper right, the same as the upper two 
panels of Fig. [9] except for the following. In the upper left panel the integrated intervals for the blue lobe 
and grey-scale image are 3.5 to 5.5 km s -1 and 6.5 to 7.5 km s _1 , respectively. The blue contour levels are 
50%, 60%,..., 90% of the peak value. In the upper right panel the position angle is 0° and contour levels 
are 0.5 to 2 K by 0.3 K, 2.5 to 6.5 K by 0.5 K. middle, the same as the lower panel of Fig. [9] lower left, 
13 CO integrated intensity map overlaid on the 13 CO grey-scale image, the integrated intervals for the blue 
lobe and the grey-scale image and the blue contour levels are the same with those in the upper left panel of 
this figure, lower right, P-V diagram of 13 CO, through the same position angle in the upper left panel of 
this figure. Contour levels are 0.15 to 0.9 K by 0.15 K, 1.1 to 1.9 K by 0.2 K and 2.1 to 3.9 K by 0.3 K. 
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Fig. 30.- TMO_26 (SST 042318.2+264115): upper left and upper right, the same as the upper two 
panels of Fig. [S] except for the following. In the upper left panel the integrated intervals for the red lobe 
and grey-scale image are 7.5 to 10 km s -1 and 6 to 7 km s _1 , respectively. The red contour levels are 40%, 
50%,..., 90% of the peak value. In the upper right panel the position angle is 50°. middle, the same as the 
lower panel of Fig. [6] lower left, 13 CO integrated intensity map overlaid on the 13 CO grey-scale image, the 
integrated intervals for the red lobe and grey-scale image are the same with those in the upper left panel of 
this figure. The red contour levels are 30%, 40%,..., 90% of the peak value, lower right, the same as the 
lower right panel of Fig. [Till except that the contour levels are 0.15 to 0.9 K by 0.15 K, 1.1 to 2.3 K by 0.2 K. 
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Fig. 31. TMO_27 (SST 042656.2+244335): The same as Fig. (30]except for the following. In the upper left 
panel the integrated intervals for the red lobe and grey-scale image are 8.3 to 9.7 km s -1 and 6 to 7.5 km 
s -1 , respectively. In the upper right panel the position angle is 40°. In the lower left panel the integrated 
intervals for the red lobe and grey-scale image and the red contour levels are the same with those in the 
upper left panel of this figure. In the lower right panel the position angle is the same with the upper right 
panel of this figure and the contour levels are 0.15 to 0.9 K by 0.15 K, 1.1 to 2.1 K by 0.2 K. 
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Fig. 32.- TMO_28 (SST 042702.6+260530): upper left and upper right, the same as the upper two 
panels of Fig. [5] except for the following. In the upper left panel the integrated intervals for the blue and red 
lobe are 3.4 to 5 km s -1 and 8 to 13 km s -1 , respectively. The blue and red contour levels are both 50%, 
60%,..., 90% of their peak value. In the upper right panel the position angle is 90° and contour levels are 0.7 
to 1.6 K by 0.3 K, 2 to 4.5 K by 0.5 K. middle left and middle right, the same as the lower two panels 
of Fig. [5] lower left, 13 CO integrated intensity map overlaid on the 13 CO grey-scale image, the integrated 
intervals for the blue lobe, red lobe and grey-scale image and the blue and red contour levels are the same 
with those in the upper left panel of this figure, lower right, the same as the lower right panel of Fig. [5H 
except for that the position angle is 90° and the contour levels are 0.2 to 2 K by 0.3 K, 2.5 to 4.5 K by 0.5 
K. 
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Fig. 33.- TMO_29 (SST 042702.8+254222): The same as Fig. [32] except for the following. In the upper 
left panel the integrated intervals for the blue lobe, red lobe and grey-scale image are 3.5 to 5.5 km s _1 , 8 to 
9.5 km s" 1 and 6 to 7 km s _1 , respectively. The blue and red contour levels are 60%, 70%,..., 90% and 30%, 
40%,..., 90% of their peak value, respectively. In the upper right panel the position angle is 160° and the 
contour levels are 0.4 to 1.6 K by 0.3 K, 2 to 4.5 K by 0.5 K. In the lower left panel the integrated interval 
for the blue lobe is 4.5 to 5.9 km s -1 . The blue levels are 40%, 50%,..., 90% of the peak value. In the lower 
right panel the position angle is 75° and the contour levels are 0.15 to 0.9 K by 0.15 K, 1.1 to 1.9 K by 0.2 
K and 2.2 to 2.8 K by 0.3 K. 
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Fig. 34.— TMO_30 (SST 042757.3+261918): The same as Fig. [T7] except for the following. In the upper- 
left panel the integrated intervals for the red lobe and grey-scale image are 8 to 10.5 km s -1 and 6 to 7.5 km 
s _1 , respectively. The red contour levels are 30%, 40%,..., 90% of the peak value. In the upper right panel 
the position angle is 90° and the contour levels are 0.4 to 1.6 K by 0.3 K, 2 to 7 K by 0.5 K. 
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Fig. 35- TM0_31 (SST 042810.4+243553): The same as Fig. [T7] except for the following. In the upper- 
left panel the integrated intervals for the red lobe and grey-scale image are 8 to 9.5 km s -1 and 6 to 7 km 
s _1 , respectively. The red contour levels are 30%, 45%,..., 90% of the peak value. In the upper right panel 
the position angle is 160° and the contour levels are 0.4 to 1.6 K by 0.3 K, 2 to 5 K by 0.5 K. 
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Fig. 36.- TMO_32 (SST 043051.7+244147): The same as Fig. [T7] except for the following. In the upper 
left panel the integrated intervals for the red lobe and grey-scale image are 8 to 11.3 km s -1 and 5.5 to 7 km 
s _1 , respectively. The red contour levels are 30%, 45%,..., 90% of the peak value. In the upper right panel 
the position angle is 60° and the contour levels are 0.4 to 1.6 K by 0.3 K, 2 to 7 K by 0.5 K. 
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Fig. 37.- TMO_33 (SST 043215.4+242859): The same as Fig. [^except for the following. In the upper left 
panel the integrated intervals for the blue lobe, red lobe and grey-scale image are 3 to 4 km s -1 , 8 to 10.5 
km s _1 and 5 to 7 km s^ 1 , respectively. In the upper right panel the position angle is 30° and the contour 
levels are 0.4 to 1.6 K by 0.3 K, 2 to 5.5 K by 0.5 K. In the lower right panel the contour levels are 0.15 to 
0.9 K by 0.15 K, 1.1 to 3.2 K by 0.2 K. 
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Fig. 38.- TMO_34 (SST 043307.8+261606): upper left and upper right, the same as the upper two 
panels of Fig. [9] except for the following. In the upper left panel the integrated intervals for the blue lobe 
and grey-scale image are 2.3 to 3.5 km s -1 and 5.5 to 7 km s _1 , respectively. The blue contour levels are 
30%, 40%,..., 90% of the peak value. In the upper right panel the position angle is 130° and contour levels 
are 0.4 to 1.6 K by 0.3 K, 2 to 4 K by 0.5 K. lower left, the same as the lower panel of Fig. [9] lower right, 
the same as the lower left panel of Fig. [25] 
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Fig. 39. TMO_35 (SST 043310.0+243343): The same as Fig. [T7] except for the following. In the upper- 
left panel the integrated intervals for the red lobe and grey-scale image are 8 to 9.5 km s -1 and 6 to 7 km 
s _1 , respectively. The red contour levels are 50%, 60%,..., 90% of the peak value. In the upper right panel 
the position angle is 140° and the contour levels are 0.4 to 1.6 K by 0.3 K, 2 to 5.5 K by 0.5 K. 
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Fig. 40.- TMO_36 (SST 043316.5+225320): The same as Fig. [32] except for the following. In the upper 
left panel the integrated intervals for the blue lobe, red lobe and grey-scale image are 2 to 4 km s -1 , 7 to 8 
km s" 1 and 4.5 to 6.5 km s -1 , respectively. In the upper right panel the position angle is 0° and the contour 
levels are 0.7 to 1.6 K by 0.3 K, 2 to 4 K by 0.5 K. In the lower right panel the contour levels are 0.3 to 1.2 
IC by 0.3 K, 1.5 to 6.5 IC by 0.5 K. 
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Fig. 41. TMO_37 (SST 043334.0+242117): The same as Fig. [T7] except for the following. In the upper- 
left panel the integrated intervals for the red lobe and grey-scale image are 8 to 9.5 km s -1 and 6 to 7 km 
s _1 , respectively. The red contour levels are 30%, 45%,..., 90% of the peak value. In the upper right panel 
the position angle is 35° and the contour levels are 0.4 to 1.6 K by 0.3 K, 2 to 5 K by 0.5 K. 
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Fig. 42.— TMO_38 (SST 043336.7+260949): The same as Fig. [38] except for the following. In the upper- 
left panel the integrated intervals for the blue lobe and grey-scale image are 3.0 to 4.1 km s -1 and 5.5 to 6.5 
km s _1 , respectively. The blue contour levels are 30%, 45%,..., 90% of the peak value. In the upper right 
panel the position angle is 70°. 
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Fig. 43.- TMO_39 (SST 043557.6+225357): The same as Fig. [32] except for the following. In the upper 
left panel the integrated intervals for the blue lobe, red lobe and grey-scale image are 2 to 4 km s _1 , 7 to 8 
km s -1 and 5 to 6 km s _1 , respectively. The blue and red contour levels are 35%, 50%,..., 95% and 60%, 
70%,..., 90% of their peak value, respectively. In the upper right panel the contour levels are 0.7 to 1.6 K 
by 0.3 K, 2 to 7 K by 0.5 K. In the lower right panel the contour levels are 0.15 to 0.9 K by 0.15 K, 1.1 to 
1.9 K by 0.2 K and 2.2 to 4 K by 0.3 K. 
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Fig. 44.— TMO_40: The same as Fig. [IT] except for the following. In the upper left panel the integrated 
intervals for the red lobe and grey-scale image are 8.5 to 11 km s -1 and 4.5 to 7 km s _1 , respectively. The 
red contour levels are 30%, 40%,..., 90% of the peak value. In the upper right panel the position angle is 90° 
and the contour levels are 0.7 to 1.6 K by 0.3 K, 2 to 5 K by 0.5 K. 
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Fig. 45.- TM0_41 (SST 043913.8+255320): The same as Fig. [32] except for the following. In the upper 

left panel the integrated intervals for the blue lobe and grey-scale image are 1.3 to 4.5 km s -1 and 5 to 6 
km s _1 , respectively. The blue and red contour levels are 45%, 60%,..., 90% and 30%, 45%,..., 90% of their 
peak value, respectively. In the upper right panel the position angle is 52° and the contour levels are 0.4 to 
1.6 K by 0.3 K, 2 to 5 K by 0.5 K. In the lower left panel the blue and red contour levels are 60%, 75%, 90% 
and 45%, 60%,..., 90% of their peak value, respectively. In the lower right panel the contour levels are 0.1 
to 1.6 K by 0.3 K, 2 to 5 K by 0.5 K. 
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Fig. 46.- TMO_42 (SST 044802.3+253359): The same as Fig. [32] except for the following. In the upper 
left panel the integrated intervals for the blue lobe, red lobe and grey-scale image are 2.8 to 3.5 km s _1 , 8 
to 9.5 km s -1 and 5 to 7 km s _1 , respectively. The blue and red contour levels are both 30%, 45%,..., 90% 
of their peak value. In the upper right panel the position angle is 100° and the contour levels are 0.4 to 1.6 
K by 0.3 K, 2 to 5 K by 0.5 K. In the lower right panel the contour levels are 0.1 to 1.6 K by 0.3 K, 2 to 6.5 
K by 0.5 K. 
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Fig. 47. TMO_43 (SST 041851.4+282026): The same as Fig. [9] except for the following. In the upper left 
panel the integrated intervals for the blue lobe and grey-scale image are 4 to 5 km s -1 and 6.5 to 8 km s _1 , 
respectively. The blue contour levels are 40%, 50%,..., 90% of the peak value. In the upper right panel the 
position angle is 135° and contour levels are 0.4 to 1.6 K by 0.3 K, 2 to 6.5 K by 0.5 K. 
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Fig. 48.— TMO_44 (SST 042021.4+281349): The same as Fig. [5]except for the following. In the upper left 
panel the integrated intervals for the blue lobe,red lobe and grey-scale image are 4.5 to 5.5 km s _1 , 8.3 to 
8.8 km s” 1 and 6.5 to 7.5 km s _1 , respectively. The blue and red contour levels are both 30%, 40%,..., 90% 
of their peak value. In the upper right panel the position angle is 165° and contour levels are 0.4 to 5.5 K 
by 0.3 K. 
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Fig. 49.— TMO_45 (SST 042653.3+255858): The same as Fig. [5] except for the following. In the upper left 
panel the integrated intervals for the blue lobe,red lobe and grey-scale image are 2 to 5.5 km s _1 , 7.5 to 10 
km s" 1 and 6 to 7 km s _1 , respectively. The blue and red contour levels are both 60%, 70%,..., 90% of their 
peak value. In the upper right panel the position angle is 60° and contour levels are 0.4 to 1.6 K by 0.3 K, 
2 to 3.5 K by 0.5 K. 
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Fig. 50. TMO_46: The same as Fig. [5] except for the following. In the upper left panel the integrated 

intervals for the blue lobe,red lobe and grey-scale image are 1 to 4 km s _1 , 7.5 to 9 km s -1 and 4.5 to 7 km 
s _1 , respectively. The blue and red contour levels are 30%, 45%,..., 90% and 60%, 70%,..., 90% of their peak 
value, respectively. In the upper right panel the position angle is 20° and contour levels are 0.4 to 1.6 K by 
0.3 K, 2 to 4.5 K by 0.5 K. 
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Fig. 51. TMO_47 (SST 043535.3+240819): The same as Fig. [35] except for the following. In the upper 
left panel the integrated intervals for the blue lobe, red lobe and grey-scale image are 4 to 4.8 km s -1 , 7 to 
8.5 km s -1 and 5 to 7 km s _1 , respectively. The red contour levels are 30%, 45%,..., 90% of the peak value. 
In the upper right panel the position angle is 0° and the contour levels are 0.7 to 1.6 K by 0.3 K, 2 to 6 K 
by 0.5 K. In the lower right panel the contour levels are 0.1 to 1.6 K by 0.3 K, 2 to 4 K by 0.5 K. 
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Fig. 52.— TMO_48 (SST 041535.6+284741): The same as Fig. [5] except for the following. In the upper left 
panel the integrated intervals for the blue lobe,red lobe and grey-scale image are -1 to 4 km s _1 , 9.25 to 13 
km s -1 and 6.5 to 8.5 km s _1 , respectively. The red contour levels are 40%, 50%,..., 90% of the peak value. 
In the upper right panel the position angle is 0° and contour levels are 0.4 to 1.6 K by 0.3 K, 2 to 6.5 K by 
0.5 K. 
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Fig. 53.- TMO_49 (SST 041733.7+282046): The same as Fig. [6] except for the following. In the upper left 
panel the integrated interval for the red lobe is 9.5 to 13 km s -1 . In the upper right panel the position angle 
is 23° and contour levels are 0.4 to 1.6 K by 0.3 K, 2 to 7 K by 0.5 K. 
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Fig. 54.— TMO_50 (SST 041810.5+284447): The same as Fig. [6] except for the following. In the upper left 
panel the integrated intervals for the red lobe and grey-scale image are 9.3 to 10 km s -1 and 6.5 to 8.5 km 
s _1 , respectively. The red contour levels are 50%, 60%,..., 90% of the peak value. In the upper right panel 
the position angle is 30° and contour levels are 0.4 to 1.6 K by 0.3 K, 2 to 5 K by 0.5 K. 
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Fig. 55.— TM0_51 (SST 041831.1+281629): The same as Fig. [6]cxcept for the following. In the upper left 
panel the integrated intervals for the red lobe and grey-scale image are 9.3 to 10.2 km s -1 and 6.5 to 8 km 
s _1 , respectively. The red contour levels are 30%, 40%,..., 90% of the peak value. In the upper right panel 
the position angle is 80° and contour levels are 0.5 to 2 K by 0.3 K, 2.5 to 6 K by 0.5 K. 


75 


Fig. 56.— TMO_52 (SST 041831.2+282617): The same as Fig. [6]cxcept for the following. In the upper left 
panel the integrated intervals for the red lobe and grey-scale image are 9.3 to 9.8 km s -1 and 6 to 8 km s _1 , 
respectively. The red contour levels are 30%, 40%,..., 90% of the peak value. In the upper right panel the 
position angle is 60° and contour levels are 0.4 to 1.6 K by 0.3 K, 2 to 6 K by 0.5 K. 
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Fig. 57. TMO_53 (SST 041841.3+282725): The same as Fig. [6] except for the following. In the upper left 
panel the integrated intervals for the red lobe and grey-scale image are 9.3 to 9.9 km s -1 and 6 to 8 km s _1 , 
respectively. The red contour levels are 60%, 70%,..., 90% of the peak value. In the upper right panel the 
position angle is 164° and contour levels are 0.4 to 1.6 K by 0.3 K, 2 to 6 K by 0.5 K. 
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Fig. 58.— TMO_54 (SST 042154.5+265231): The same as Fig. [5] except for the following. In the upper left 
panel the integrated intervals for the blue lobe,red lobe and grey-scale image are 3.8 to 4.5 km s _1 , 7.5 to 
8.3 km s -1 and 5.5 to 7 km s _1 , respectively. The blue and red contour levels are 30%, 45%,..., 90% and 
50%, 60%,..., 90% of their peak value, respectively. In the upper right panel the position angle is 0° and 
contour levels are 0.5 to 2 K by 0.3 K, 2.5 to 4 K by 0.5 K. 
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Fig. 59.— TMO_55 (SST 042904.9+264907): The same as Fig. [6] except for the following. In the upper left 
panel the integrated intervals for the red lobe and grey-scale image are 8.3 to 9 km s^ 1 and 5.5 to 7 km s -1 , 
respectively. In the upper right panel the position angle is 50° and contour levels are 0.4 to 4.2 K by 0.3 K. 
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Fig. 60.- TMB_01: upper left, 13 CO integrated intensity map, integrated over 7.3 to 8.4 km s _1 . The 
starting contour and the contour step are 3.3 K km s^ 1 and 0.6 K km s -1 , respectively. The black dashed 
line represents a cut for the P-V diagram shown in the upper middle panel, and the amaranth open circle 
indicates the origin on the Y-axis of the upper middle panel. The purple dashed circle shows the approximate 
extent of the CO emission of the bubble. The black open circle on the lower left corner marks the beam, 
upper middle, P-V diagram of 13 CO, through the slice shown by the black dashed line in the upper left 
panel at a position angle of 110°. The contour levels are 0.1 to 1.2 K by 0.1 K. The two black dashed vertical 
lines show the expansion velocity interval determined by visual inspection, upper right, gaussian fit to the 
azimuthally averaged profile of 13 CO intensity, lower, Channel maps of 13 CO emission, the number on the 
upper right corner of each panel indicates the central velocity of the channel map. The purple dashed circles 
are the same as that in the upper left panel. 
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Fig. 61.— TMB_02: The same as Fig. [60] except for the following. In the upper left panel the integrated 
velocity interval is from 5.4 to 6.5 km s -1 , the starting contour and the contour step are 2.7 K km s -1 and 
0.5 K km s~\ respectively. In the upper middle panel the position angle is 80° and the contour levels are 
0.1 to 1.0 K by 0.1 K. 


Fig. 62.— TMB_03: The same as Fig. [60] except for the following. In the upper left panel the integrated 
velocity interval is from 4.6 to 6.2 km s -1 and the contour step is 0.9 K km s” 1 . In the upper middle panel 
the position angle is 160° and the contour levels are 0.1 to 1.7 K by 0.2 K. 


81 


Fig. 63.- TMB_04: The same as Fig. [60] except for the following. In the upper left panel the integrated 
velocity interval is from 6.0 to 7.0 km s -1 and the contour step is 0.8 K km s -1 . In the upper middle panel 
the position angle is 20° and the contour levels are 0.1 to 2.2 K by 0.2 K. 

Fig. 64.— TMB_05: The same as Fig. 1501 except for the following. In the upper left panel the integrated 
velocity interval is from 4.4 to 6.0 km s -1 , the starting contour and the contour step is 2.7 K km s -1 and 1 
K km s -1 , respectively. In the upper middle panel the position angle is 10° and the contour levels are 0.1 
to 1.9 K by 0.2 K. 

Fig. 65.— TMB_06: The same as Fig. [Ml except for the following. In the upper left panel the integrated 
velocity interval is from 5.2 to 6.2 km s -1 , the starting contour and the contour step are 2.1 K km s -1 and 
0.3 K km s -1 , respectively. 

Fig. 66.- TMB_07: The same as Fig. [Ml except for the following. In the upper left panel the integrated 
velocity interval is from 5.7 to 7.8 km s -1 , the starting contour and the contour step are 4.7 K km s -1 and 
0.9 K km s -1 , respectively. In the upper middle panel the position angle is 10° and the contour levels are 
0.1 to 1.1 K by 0.1 K. 

Fig. 67. TMB_08: The same as Fig. [Ml except for the following. In the upper left panel the integrated 

velocity interval is from 5.7 to 7.8 km s^ 1 and the contour step is 1 K km s -1 . In the upper middle panel 
the position angle is 10° and the contour levels are 0.1 to 1.1 K by 0.1 K. 

Fig. 68.— TMB_09: The same as Fig. [Ml except for the following. In the upper left panel the integrated 
velocity interval is from 6.2 to 7.0 km s 1 , the starting contour and the contour step are 3.4 K km s -1 and 
0.7 K km s -1 , respectively. In the upper middle panel the position angle is 80°. 

Fig. 69.- TMB_10: The same as Fig. [Ml except for the following. In the upper left panel the integrated 
velocity interval is from 4.4 to 5.7 km s -1 , the starting contour and the contour step are 2.3 K km s -1 and 
0.8 K km s -1 , respectively. In the upper middle panel the contour levels are 0.1 to 1.8 K by 0.2 K. 

Fig. 70. TMB_11: The same as Fig. [Ml except for the following. In the upper left panel the integrated 

velocity interval is from 6.2 to 8.1 km s -1 , the starting contour and the contour step are 4.8 K km s -1 and 
0.9 K km s -1 , respectively. In the upper middle panel the position angle is 40° and the contour levels are 
0.1 to 1.6 K by 0.2 K. 

Fig. 71.— TMB_12: The same as Fig. [Ml except for the following. In the upper left panel the integrated 
velocity interval is from 5.4 to 7.8 km s _1 , the starting contour and the contour step are 6.2 K km s -1 and 
0.9 K km s _1 , respectively. In the upper middle panel the position angle is 160° and the contour levels are 
0.1 to 1.5 K by 0.2 K. 

Fig. 72.— TMB_13: The same as Fig. [Ml except for the following. In the upper left panel the integrated 
velocity interval is from 4.4 to 6.5 km s — , the starting contour and the contour step are 5.0 K km s -1 and 
0.9 K km s -1 , respectively. In the upper middle panel the position angle is 40° and the contour levels are 
0.1 to 1.7 K by 0.2 K. 

Fig. 73.- TMB_14: The same as Fig. [Ml except for the following. In the upper left panel the integrated 
velocity interval is from 0.7 to 2.5 km s” 1 and the starting contour is 1.3 K km s -1 . In the upper middle 
panel the position angle is 130° and the contour levels are 0.1 to 0.8 K by 0.1 K. 
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Fig. 74.— TMB_15: The same as Fig. [60] except for the following. In the upper left panel the integrated 
velocity interval is from 6.8 to 7.8 km s -1 , the starting contour and the contour step are 2.9 K km s -1 and 
0.8 K km s -1 , respectively. In the upper middle panel the position angle is 20° and the contour levels are 
0.1 to 1.4 K by 0.1 K. 


Fig. 75.— TMB_16: The same as Fig. [60] except for the following. In the upper left panel the integrated 
velocity interval is from 6.5 to 8.1 km s -1 , the starting contour and the contour step are 4.2 K km s -1 and 
0.9 K km s -1 , respectively. In the upper middle panel the position angle is 10° and the contour levels are 
0.1 to 1.8 K by 0.2 K. 


Fig. 76.- TMB_17: The same as Fig. [60] except for the following. In the upper left panel the integrated 
velocity interval is from 5.2 to 6.0 km s' 1 and the starting contour is 1.6 K km s . In the upper middle 
panel the position angle is 170° and the contour levels are 0.1 to 1.3 K by 0.1 K. 


Fig. 77. TMB_18: The same as Fig. [60] except for the following. In the upper left panel the integrated 

velocity interval is from 5.2 to 7.0 km s _1 , the starting contour and the contour step are 7.2 K km s -1 and 
1 K km s' 1 , respectively. In the upper middle panel the position angle is 170° and the contour levels are 
0.1 to 1.9 K by 0.2 K. 


Fig. 78.- TMB_19: The same as Fig. [60] except for the following. In the upper left panel the integrated 
velocity interval is from 4.9 to 6.8 km s _1 , the starting contour and the contour step are 5.0 K km s _1 and 
0.8 K km s _1 , respectively. In the upper middle panel the contour levels are 0.1 to 1.5 K by 0.1 K. 


Fig. 79.— TMB_20: The same as Fig. [001 except for the following. In the upper left panel the integrated 
velocity interval is from 5.4 to 6.5 km s' 1 and the starting contour is 3.6 K km s . In the upper middle 
panel the position angle is 130° and the contour levels are 0.1 to 1.8 K by 0.2 K. 


83 


Fig. 80. TMB_21: The same as Fig. ISU1 except for the following. In the upper left panel the integrated 

velocity interval is from 1.4 to 3.6 km s' 1 , the starting contour and the contour step are 2.2 K km s' 1 and 
0.7 K km s' 1 , respectively. In the upper middle panel the position angle is 140° and the contour levels are 
0.1 to 0.8 K by 0.1 K. 


Fig. 81.— - TMB_22: The same as Fig. 1501 except for the following. In the upper left panel the integrated 
velocity interval is from 4.4 to 5.4 km s' 1 , the starting contour and the contour step are 2.0 K km s' 1 and 
0.7 K km s' 1 , respectively. In the upper middle panel the position angle is 120° and the contour levels are 
0.1 to 1.5 K by 0.1 K. 


Fig. 82.— TMB_23: The same as Fig. ISU1 except for the following. In the upper left panel the integrated 
velocity interval is from 4.9 to 6.0 km s' 1 , the starting contour and the contour step are 3.6 K km s' 1 and 
0.4 K km s' 1 , respectively. In the upper middle panel the position angle is 120° and the contour levels are 
0.1 to 1.7 K by 0.2 K. 


Fig. 83.— TMB_24: The same as Fig. [60] except for the following. In the upper left panel the integrated 
velocity interval is from 4.9 to 6.0 km s' 1 , the starting contour and the contour step are 4.2 K km s' 1 and 
0.5 K km s' 1 , respectively. In the upper middle panel the position angle is 70° and the contour levels are 
0.1 to 0.9 K by 0.1 K. 


Fig. 84.— TMB_25: The same as Fig. [60] except for the following. In the upper left panel the integrated 
velocity interval is from 4.9 to 7.8 km s' 1 , the starting contour and the contour step are 8.2 K km s' 1 and 
1.2 K km s' 1 , respectively. In the upper middle panel the position angle is 150° and the contour levels are 
0.1 to 1.0 K by 0.1 K. 


Fig. 85.~ TMB_26: The same as Fig. [SOI except for the following. In the upper left panel the integrated 
velocity interval is from 2.0 to 4.1 km s' 1 , the starting contour and the contour step are 1.5 K km s' 1 and 
0.8 K km s' 1 , respectively. In the upper middle panel the position angle is 80° and the contour levels are 
0.1 to 0.6K by 0.1 K. 


Fig. 86.- TMB_27: The same as Fig. [60] except for the following. In the upper left panel the integrated 
velocity interval is from 2.5 to 4.9 km s' 1 and the starting contour is 1.8 K km s' 1 . In the upper middle 
panel the position angle is 140°. 


Fig. 87. TMB_28: The same as Fig. [60] except for the following. In the upper left panel the integrated 

velocity interval is from 4.4 to 7.6 km s' 1 , the starting contour and the contour step are 7.5 K km s' 1 and 
1.2 K km s' 1 , respectively. In the upper middle panel the position angle is 40°. 
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Fig. 88.— TMB_29: The same as Fig. [60] except for the following. In the upper left panel the integrated 
velocity interval is from 4.4 to 5.7 km s -1 , the starting contour and the contour step are 4.0 K km s -1 and 
0.4 K km s -1 , respectively. In the upper middle panel the position angle is 20° and the contour levels are 
0.1 to 1.3 K by 0.1 K. 


Fig. 89.- TMB_30: The same as Fig. 1501 except for the following. In the upper left panel the integrated 
velocity interval is from 4.4 to 6.0 km s _1 and the starting contour is 5.1 K km s -1 . In the upper middle 
panel the position angle is 130° and the contour levels are 0.1 to 1.4 K by 0.1 K. 


Fig. 90. TMB_31: The same as Fig. [60] except for the following. In the upper left panel the integrated 

velocity interval is from 4.6 to 6.0 km s _1 and the starting contour is 1.7 K km s . In the upper middle 
panel the position angle is 120° and the contour levels are 0.1 to 0.8 K by 0.1 K. 


Fig. 91.— - TMB_32: The same as Fig. [60] except for the following. In the upper left panel the integrated 
velocity interval is from 2.2 to 4.1 km s -1 , the starting contour and the contour step are 1.4 K km s -1 and 
0.4 K km s -1 , respectively. In the upper middle panel the position angle is 50° and the contour levels are 
0.1 to 0.5 K by 0.1 K. 


Fig. 92.— TMB_33: The same as Fig. [60] except for the following. In the upper left panel the integrated 
velocity interval is from 7.0 to 7.8 km s -1 , the starting contour and the contour step are 2.0 K km s -1 and 
0.4 K km s“\ respectively. In the upper middle panel the position angle is 50° and the contour levels are 
0.1 to 2.2 K by 0.2 K. 


Fig. 93.— TMB_34: The same as Fig. 1501 except for the following. In the upper left panel the integrated 
velocity interval is from 5.2 to 7.8 km s 1 , the starting contour and the contour step are 3.8 K km s -1 and 
0.7 K km s -1 , respectively. In the upper middle panel the position angle is 40°. 


Fig. 94.— TMS_35: The same as Fig. IMl except for the following. In the upper left panel the integrated 
velocity interval is from 4.9 to 6.2 km s _1 and the starting contour is 4.0 K km s -1 . In the upper middle 
panel the position angle is 50° and the contour levels are 0.1 to 0.9 K by 0.1 K. 


Fig. 95.— TMB_36: The same as Fig. [SOI except for the following. In the upper left panel the integrated 
velocity interval is from 5.7 to 7.3 km s -1 , the starting contour and the contour step are 4.0 K km s -1 and 
0.4 K km s -1 , respectively. In the upper middle panel the position angle is 30° and the contour levels are 
0.1 to 0.8 K by 0.1 K. 


Fig. 96.- TMB_37: The same as Fig. [60] except for the following. In the upper left panel the integrated 
velocity interval is from 3.3 to 5.4 km s -1 , the starting contour and the contour step are 2.2 K km s -1 and 
0.7 K km s -1 , respectively. In the upper middle panel the contour levels are 0.1 to 1.1 K by 0.1 K. 
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A. Derivations of Outflow Parameters 

To calculate molecular outflow parameters we need first to obtain the column density first. A simple 
solution of the equation of radiative transfer is 


T s = [T ex -T bg ]( l-e-^), 


(Al) 


where T s is the source temperature, T ex is the excitation temperature and T bg is the background temperature, 
and v is the frequency of the transition. The modified Planck function J is defined as 


J(T) = 


hv/k 

hZ 7 1 

e -pr — 1 


(A2) 


where k is Boltzmann’s constant and h is Pl anck’s constant. The definition of optical depth in terms of 
upper level column density is expressed as (Wilson et ah 201 3) 




(A3) 


Assuming r <C 1, hv/k -C T ex and T bg <C T ex , we get the column density of the rotational upper level of the 
transition in the outflow by combining Eq. m with Eq. 


iV u ( 12 CO) = 


87 xkv 2 
hc 3 A u \ 


T R dv, 


(A4) 


where T s is the observed source antenna temperature with proper correction for antenna efficiency, c is the 
speed of light, and A u i is the spontaneous transition rate from the upper level (J + 1 ) to the lower level(J), 
which can be expressed as 

64t J+l 


Aul 3 he 3 2 J + 3 ’ 

where gd is the permanent electric dipole moment of a molecule and J = 0. 
The total column density of the outflow is 


(AS) 


NtotC 2 CO) = 


iV u ( 12 CO) 

/u ’ 


(A 6 ) 


where / u is the fraction of the 12 CO in the upper level of the transition. Under local thermal equilibrium(LTE) 
conditions, f u is given by 

g u exp(-hi//kT ex ) 


/» = 


Q(Tex) 


(A7) 


where the statistical weight of the upper level g u = 2 (J +1) +1, the LTE partition function (for k T PX h B ) 
Q(T e x ) = kT ex /hB and the rotational constant B = v/\2 (J + 1)] for the J + 1 —>• J transition ([Tennyson 
20051) . Then we can derive the total column density of outflow from Eq. (IA4I) . Eq. (IASI) . Eq. (IA6D and Eq. 


C3, 


NtotC z CO) = 


47 T 3 ii\hv 2 exp(— hv/kT ex ) 


T s dv 


(A 8 ) 


If there is a high velocity wing in 12 CO but not in 13 CO profile, we assume 12 CO is optically thin. Then we 
can calculate the column density of outflow from Eq. (1A8[) . If there is high velocity wing both in 12 CO and 
13 CO profile, we can correct the optical depth of 12 CO using the following equation 


T*( 12 CO) l-e- r ( 12c °) 

T*( 13 CO) “ l - e -T( 13 c°) • 


(A9) 
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Here T*( 12 CO) and T*( 13 CO) are the antenna temperatures of 12 C0 and 13 C0 (with proper correction 
for antenna efficiency), respectively. r( 12 CO) and r( 13 CO) are the optical dept hs of 12 C0 and 13 CO, 
respectively. We assume the abundance ratio of 12 C0 to 13 CO is 65 ( Langer fc Penziaslll993l f. The correction 
factor for opacity is defined as 

f - T{12 ° 0) (A10) 

J' T ~ l _ g-r( 12 CO) ' 

Then we get the corrected total column density of outflow as 


iVctot( 12 CO) = f T N t ot ( 12 CO). 


(All) 


After obtaining the column density we can calculate other parameters of outflow. The mass of outflow can 
be calculated from 


M gas = Wot( 12 CO)[H 2 /COKm(H)S, 


(A12) 


where /i g =2.72 is the mean molecular weight ( Bruntll2010lf . ?n(H)=1.67 x 10~ 24 g is the mass of a hydrogen 
atom, and [H 2 ]/[CO] is assumed to be 10 4 , and S is the area of the outflow. 

The momentum (P) and energy ( E ) of the outflow can be calculated from 


P = M gas \v\, (A13) 

E=-M Eas v\ (A14) 

where v is the average velocity of the outflow relative to the cloud systemic velocity and M gas is obtained 
from Ea. (IA12|) . 

The dynamical timescale fdyn can be estimated from 


td 


yn — 



(A15) 


where L is the typical linear scale of the outflow lobe. The outflow luminosity, Pfl OW ) can be estimated by 
dividing the kinetic energy by the dynamical timescale. It can be expressed as 
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